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Chapter 1 



Introduction 



We begin with a description of recent numerical and analytical results that are 
closely related to the results of this paper. 



In 1978 Holian and Straub [HS| conducted an extensive series of numerical ex- 



periments on a driven, semi-infinite lattice 

(1.1) Xn = F{Xn~l - Xn) - F{Xn - Xn+l), n = 1, 2, . . . , 

with initial conditions 

(1.2) Xn{0) = nd, i;„(0) = 0, n = l,2,..., d constant, 

for a variety of force laws F, and in the case that the velocity of the driving particle 
xo is fixed, 

(1.3) xo{t) = 2at, t>0, a > 0. 

They discovered, in particular, a striking new phenomenon - the existence of a 
critical "shock" strength fflcrit- If a < acrit = acrit(^)) then in the frame moving 



with the particle xq, they observed behavior similar to that shown in Figure |1.4 . 
Thus the particles come to rest in a regular lattice behind the driver. However, if 
> Ocrit = cLcrit{F), then, again in the frame of the driver, they observed behavior as 



in Figure 1.5. Now the particles do not come to rest behind the driver, but execute 
an on-going binary oscillation (i.e. Xn{t + T) = Xn{t), Xn{t) = x„+2(t) + const.). 
This is a marvelous, fundamentally nonlinear phenomenon; if F is linear, the effect 
is absent. 

This phenomenon has now been observed for many different singular and non- 
singular, nonlinear force laws F, but an explanation of the phenomenon from first 
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Figure 1.4: Motion of the first ten particles of a lattice described by the above system 
(1-1) - (1-3) with F(x) = e^,c/ = 0,a = .5, in the frame of xq (case a < Ocrit)- 
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Figure 1.5: Motion of the first ten particles of a lattice described by the above system 
(1-1) - (1-3) with F(x) = e^,d = 0,a = 2, in the frame of xq (case a > acritj- 
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principles in the general case has not yet been given. We believe that the phe- 
nomenon is present for a very wide class of genuinely nonlinear F (in particular, if 
F" > 0) with F' > 0. Observe that if F' > 0, then the force on the n}^ particle 
Xn = F{xn-i — Xn) — F{xn — Xn+i) is negative for Xn > \{xn-i + Xn+i) and positive 
for Xn < \ {xn-i + Xn+i)- Thus the only equilibrium configuration is the regular 
lattice Xn = cn, c E R, and moreover, in this case, all the forces are restoring. 

In 1981, Holian, Flaschka and McLaughlin [HFM] considered the shock problem 
in the special case in which F is an exponential F{x) = e^, the so-called Toda shock 
problem. They considered this case because the Toda equation 

(1.6) Xn = e^"-i-^" - e^"-^"+i, 

with appropriate boundary conditions, is well-known to be completely integrable (a 
fact discovered by Flaschka and Manakov ||Mai] ]; see also Q) and hence there 



was the possibility of solving (IT) - (|1.3| ) explicitly and so explaining the phenomena 
observed by Holian and Straub in the special case where F{x) = e^. However, the 
driven system ( |1.1| ) - (1.3) is non-autonomous and it was not clear a priori that the 



(formal) integrability of the Toda equation could be used to analyze the system. 
For example, the one-dimensional oscillator x + ax -|- bx^ = is certainly integrable; 
however, the driven oscillator x + ax + bx^ = f{t), the so-called Duffing system, is 
far from "integrable" and requires highly sophisticated techniques for its analysis. 
Nevertheless, Holian, Flaschka and McLaughlin ([HFM|) realized that if they went 



into the frame of the driver, so that (L2), (L3) become 

(1.7) x„(0) = nd, x„(0) = -2a, n > 1, 

(T8) xo(t) = 0, 

and doubled up the system 

(1.9) Xn{t) = -x-n(.t), n < 0, 

then the full system {xn}^=-oo solves the autonomous Toda equations 

(1.10) x„ = e^"-^"^" - e^"~^"+\ -cxD < n < oo, 
with initial conditions 

(1.11) Xn{0) = dn, Xn(0) = — 2a( sgn n), — oo < n < oo. 
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But the solutions of these equations lie in a class to which the method of inverse 
scattering applies. To see what is involved we use Flaschka's variables, 



(1.12) 
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2 (-^n Xn-\-l) 



-DO < n < oo, 



and arrange these variables into a doubly-infinite tridiagonal matrix 



(1.13) 



L 



b-2 a_i 6_i O 

b-i ao bo 

bo ai bi 

O bi ••. 



which represents the state of the system at any given time, with companion matrix 



(1.14) 



B 



-b-2 6_i O 

bo 

-bo bi 

O -bi 



V 



Then, remarkably, ( 1.10| ) is equivalent to the so-called Lax-pair system 



(1.15) 



dL 
'dt 



[B,L\ = BL- LB. 



Thus the Toda equations are equivalent to an iso-spectral deformation of the matrix 
operator L. Inverse scattering theory tells us that one can solve ( 1.15| ), and hence 
( |1.1| ) - (1.3), through the scattering map for L. Rescaling time, one sees that it is 
sufficient to consider the case where the initial spacing d = {). Then at t = 0, 



(1.16) 



a sgn n, 



and one sees that the essential spectrum of L is given by two bands (cf Figure 1.17). 
The bands overlap if and only if a < 1. Holian, Flaschka and McLaughlin observed 
that supercritical behavior occurred for the Toda shock problem only if the gap was 
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Figure 1.17: The spectrum of L{0) 
open. Hence they identified Ocrit 

{F = e^) = 1. Using the inverse method they were 
able to calculate a number of other features of the Toda shock problem, such as the 
speed and the form of the shock front, and also the form of the binary oscillations. 
The problem of how to extract detailed information about the long-time behavior 
of the Toda shock problem from knowledge of the initial data using the rather 
formidable formulae of inverse scattering theory, however, remained open. 

In the early 80's, a very important development took place in the analysis of 
infinite-dimensional integrable systems in the form of the calculation by Lax and 
Levermore ([ |LL| ) of the leading order asymptotics for the zero-dispersion limit of the 
Kortweg de Vries equation, in which the weak limit of the solution as the dispersion 
coefficient tends to zero is derived and the small scale oscillations that arise are 
averaged out. This was followed in the late 80's by the calculation of Venakides 
for the higher order terms in the Lax-Levermore theory which produces the detailed 
structure of the small scale oscillations. These developments raised the possibility 
of being able to analyze the inverse scattering formulae for the solution of the Toda 



shock problem effectively as t ^ oo, and in [ VDO |, Venakides, Deift and Oba proved 
the following result in the supercritial case a > 1: 

In addition to the shock speed Vg calculated by Holian, Flaschka and McLaughlin, 
there is a second speed < vq < Vg. 

In the frame moving with the driver, as t — > oo. 



for < n/t < vq, the lattice converges to a binary oscillation Xn{t + T) 



Xnit), Xn+2it) = Xn{t)+ Constant, (cf Figure 1.5). The band structure corre- 
sponding to the asymptotic solution is [—a — 1, —a + 1] U [a — 1, a + 1]. The 
binary oscillation is connected to the driver Xo{t) = 0, through a boundary 
layer, in which the local disturbance due to the driver decays exponentially in 
n. 

ioT Vq < n/t < Vs, the asymptotic motion is a modulated, single-phase, quasi- 
periodic Toda wave with band structure [— a— 1, ^{n/t)] U [a — 1, a+1], where 
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'y{n/t) varies monotonically from — a + 1 to —a — 1 as n/t increases from vq to 



for n/t > Vs, the deviation of the particles from their initial motion —2at is 
exponentially small. The influence from the shock has not yet been felt. As 



noted in | HFM ] , for n/t ~ Vs, the motion of the lattice is described by a Toda 



solution with associated spectrum {—a — 1} U [a — 1, a + 1]. 



In 1991, again using the techniques in |LL| and Q, Kamvissis ([ |Kam| ) showed 



that in the subcritical case a < 1, in the frame moving with the driver, as t — > cxd, 
the oscillatory motion behind the shock front dies down to a quiescent regular lattice 



with spacing Xn+i — Xn ^ — 21og(l + a), (cf Figure |l.4D . 



A "Thermodynamic" Remark. 

It is easy to see that the average spacing of the binary oscillation of the asymp- 
totic state in the case a > 1 is given by — In 4a. Thus the average spacing of the 
asymptotic states is given by 

-21n(l + a), for a < 1, 
— In 4a, for a > 1. 

Observe that these expressions and their first derivatives agree at a = 1. Thus we 
may say that the density of the asymptotic state has a second order phase transition 
at a = 1. 

As in Q and 0, the ab ove results are not fully rigorous and rely on certain 
(reasonable) asymptions that have not yet been justified from first principles. In 
particular, as in [[LL[| , the contribution of the refiection coefficient associated with 
the band [a — 1, a + 1] is ignored. Also, as in 0], an ansatz is needed to control 
the long-time behavior of certain integrals. Recently in DMV| , the authors, again 



using the approach of [LL] and circumvented the first difficulty by considering 



finite dimensional approximations to the lattice of length i[t) ^ t, but they still 
need the above mentioned ansatz in order to re-derive the results in | VDQf |. In |BK|, 



Bloch and Kodama consider the Toda shock problem, both in the subcritical and 
the supercritical cases, from the point of view of Whitham modulation theory in 
which the validity of a modulated wave form for the solution is assumed a priori, 
and the parameters of the modulated wave form are calculated explicitly. More 
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recently in |]GN| , Greenberg and Nachman have considered the shock problem for a 
general force law in the weak shock limit; they are able to describe many aspects 
of the solution, including the modulated wave region where they use a KdV-type 
continuum limit. 

In a slightly different direction, motivated by the so-called von-Neumann problem 
arising in the computation of shock fronts using discrete approximations, Goodman 



and Lax |GL| and Hou and Lax [HL| observed and analyzed features strikingly 
similar to those in pS| . Finally, in an interesting series of papers spanning the 
80's, Kaup and his collaborators, [|Kau|| , pCNf , [|WK|| , use various integrable features 



of the non-autonomous system ( |1 . 1[ ) - (^^) to gain valuable insight into the Toda 
shock problem. We will return to these papers below. 



In this paper we consider the driven lattice (1.1), (1-2) in the case where the 



uniform motion of the driving particle xq is periodically perturbedlll 

{Xn = F{Xn-l - Xn) - F{Xn - Xn+l), n > 1, 
Xn{0) = ±n{0) =0, n > 1, 

xo{t) = 2at + h{jt). 

Here h{-) is periodic with period 27r and the frequency 7 > is constant. We restrict 
our attention to the case where the average value of the velocity of the driver xq = 2a 
is subcritical, i.e. a < Ccrit- (For some discussion of the supercritical case a > 1, 
see Problem 3 at the end of the Introduction below). Again we consider a variety 
of force laws F, but henceforth we restrict our attention to forces which are real 
analytic and monotone increasing in the region of interest. 
Typically we observed the following phenomena: 

In the frame moving with the average velocity 2a of the driver, as t — > 00, the 
asymptotic motion of the particles behind the shock front, is ^-periodic in time, 

(1.19) Xn{t+—)=Xnit), 0<n<^t. 

7 

Moreover, there is a sequence of thresholds, 

71 = 'yiia,h,F) > 72 = 72(a,/i,F) > • • • > 7^ = -fk{a,h,F) > • • • > 0, 
^"'"■^'^■^ 7fc — > as /c — > 00. 



^The more general initial value Xn{Q) = dn, in(0) — 0, can clearly be converted to (1.18) by 
shifting the argument of F, F{-) —> F{- ~d) : in the case of Toda, as noted above, this shift converts 
into a rescaling of the time. 
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If 7 > 7i) there exist constants c, d such that Xn—cn—d converges exponentially 
to zero as n ^ oo, (cf Figure p.6D . In other words, the effect of the oscillatory 
component of the driver does not propagate into the lattice and away from 
the boundary n = 0. The lattice behaves in a similar way to the subcritical 
case of constant driving considered by Holian, Flaschka and McLaughlin. 

If 7i > 7 > 72, then the asymptotic motion is described by a travelling wave 
(1.21) = cin + Xi(/3in + 7t), K n < t. 



transporting energy away from the driver xq. (See Figure p.7| ). Here ci = 
ci(a, /i, F, 7) and Xi{-) = Xi{- ;a,/i, F, 7) is a 27r-periodic function. 

More generally, if 7^ > 7 > 7^+1, a multi-phase wave emerges which is well- 
described by the wave form 

(1.22) Xn = CkU + XkiPin + 'yt, I32n + 2'yt,... ,l3kn + k'yt), 1 <C n <C t, 
again transporting energy away from the driver (see Figure |C.8 for the case 



k = 2). Here Ck = Ck{a, h, F,'y) and Xfc(-, . . . , •; a, /i, F, 7) is 27r-periodic in 
each of its k variables. 

Thus, at the phenomenological level, we see that the periodically driven lattice 
behaves like a long, heavy rope which one shakes up and down at one end. 



Remark: 

We have restricted our experiments to the asymptotic region 1 <^ n <^ t. However, 
we expect that the solution also exhibits many interesting phenomena when studied 
as a function of n/t. For example (see discussion on page 5), we expect that for 
7fc > 7 > 7fc+i) there will be a sequence of 2k speeds si > S2 > • • • > S2k, with the 
property that for t large, 

• for S2 < n/t < si, the solution is a modulated one-phase wave, 

• for S3 < n/t < S2, the solution is a pure one-phase wave, 

• for S4 < n/t < S3, the solution is a modulated two-phase wave, 
and so on, until 
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for S2k < n/t < S2fc-i, the solution is a modulated A;-phase wave, 



and 



for 1/t <C n/t < S2k, the region studied in this paper, we have a pure fc-phase 
wave. 



Note from the Figures C.9 - |C.ll that the above phenomena are present for both 
small and large values of the amplitude of the periodic component h of the driver. 
In the linear case, F{x) = ax, a > 0, the origin of the thresholds 71 > 72 > ... is 
simple to understand. The solution of the lattice equations 



Xn = F {Xn-1 - Xn) - F{Xn - Xn+l) = 0:{Xn+l + Xn-1 - '^Xn), n>l, 

(1 23) ^ ^ ^' " - 1' 

JTIGZ 

is easy to evaluate using Fourier methods and one sees that as t ^ 00, 
(1.24) Xn{t) — > 2a{t -n) + Y. ^"^ e^^""*' < n < 



where z^, \zra\ < 1, is the root of 



(1.25) 



{'ymY 



a 



2 z„ + 1 = 0, 



chosen, in the case \zm\ = 1, such that the energy is transported away from the 
driver. Observe that if mo > is the largest integer m for which {'ym)'^ /a — 2 < 2, 
then \zm\ = 1 for < \m\ < mo, and \zm\ < 1 for |m| > mo. Inserting this 
information into ( 1.24| ), we find that, away from the driver, an mo-phase wave 
propagates through the lattice in the region < n <^ t. Thus the threshold values 
of 7fc are given, in this case, by 



(1.26) 



Ik 



k ' 



A: = l,2,... . 



As we will see below, the above calculations are useful in understanding the asymp- 



totic state of the solution of (1.18) as t ^ 00 in the case that h is small. 

In the case of the Toda lattice, when the driving is constant the doubling-up trick 



converts the shock problem into an iso-spectral deformation (1.15) for the operator 
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L. When h is non-zero, it is no longer clear how to convert the shock problem (1.18) 
into an integrable form (although recent results of Fokas and Its |F]] suggest that 
this may still be possible to do). As a tool for analyzing ( p.. 18 ) in the Toda case we 
consider, rather, the Lax pair of operators 



(1.27) 



L 



ai 
bi 



bi 
a2 
b2 



\ 



B 



/ 6i 



V 



\ 



for the semi-infinite lattice a„ = — Xn/2, bn = \ ez^^" ^n+i) ^ n > 1. 
A straightforward calculation shows that L solves the equation 



(1.28) 



L = [B,L]-2bl{t)P, bo 



1 1 

- £2 



(xo-xi) 



which we think of as a forced Lax system. Here P = {Pij)ij>i, is a matrix operator 
with Pij = unless i = j = 1, and Pu = 1. The equation describes a motion 
that is almost, but not quite, an iso-spectral deformation of L. As t evolves, the 
essential spectrum of L{t) remains fixed, aess{L{t)) = (T(,ss{L{^)) = [a — l,a + 1], 
but eigenvalues may "leak out" from the continuum. This is true, in particular, in 
the case of constant driving xq = 2at, as was first observed by Kaup and Neuberger 



|Kg. 



In the case of constant driving with a < 1, what happens to a{L{t))l We see in 



Figure 1.29 that the eigenvalues emerge from the band [a — 1, a + 1] and eventually 
fill the larger band [—a — 1, a + 1] = [—a — 1, —a + 1] U [a — 1, a + 1]. (In the case 
a > 1, the bands [—a — 1, —a + 1], [a — 1, a + 1] are disjoint and the spectrum of L{t) 
fills these two bands separated by a gap). Thus this "ghost" band, which appeared 
as an artifact of the solution procedure through the introduction of the doubled-up 
operator L, now emerges in real form, populated by eigenvalues emerging from the 
original band [a — l,a + 1]. We learn from Figure 1.29| that for a < 1 there is no 
gap in the spectrum at t = cxd, and (hence) there are no oscillations. 

In the periodically driven case, xq = 2at + h{'yt), where 7 > 71 (and a < 1), we 
find a similar picture to Figure 1.2£ for the evolution of a{L(t)) which is displayed in 



Figure 1.30 at some later time, so that more eigenvalues are present than in Figure 
1.2£. As t ^ 00, a{L{t)) again converges to a single band and no travelling wave 



emerges. However, if 71 > 7 > 72, we find different behavior (Figure 1.31). We see 
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Figure 1.30: Evolution of 

XQ{t) = t + 0.2(sm7t + 0.5cos27i),7 = 3.1 > 71 



aim); 



driver: 
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Figure 1.31: Evolution of a{L{t)); driver: 

xo{t) = t + 0.2(sm7t + 0.5cos27t),7 = 1.8,71 > 7 > 72 

that a{L{t)) converges to two bands separated by a gap, and a single phase wave 
emerges. For 72 > 7 > 73, we see from Figure 1.32| that a{L{t)) converges to three 



bands separated by two gaps, and a two phase wave emerges, etc. 



Remarks on the eigenvalues in the gap. 

(1) The eigenvalues which can be observed in the gaps of the spectrum of the 



semi-infinite Lax operator L{t) (compare with Figures 1.31 and 1.32) are of 
two different types. They are either constant in time or they move down from 
the lower edge of one band to the upper edge of the next band. Eigenvalues 
of the second kind can be understood from the corresponding (7-gap solution 
They are connected to the zeros of a theta function, which is used in the 
construction of the g'-gap solution (see Chapter ^). Eigenvalues which are 
constant in time can be interpreted as follows. Numerical computations show 
that they correspond to eigenvectors which are moving out as t ^ cxd. Hence 
the eigenvalues do not survive in the spectrum of the limiting operator L, which 
corresponds to a lattice where all particles perform time periodic motion. In 
other words, from the spectral theoretic point of view, this is an example of the 
general phenomenon that under strong convergence of operators the spectrum 
is not necessarily conserved. 



(2) Figures 1.31 and 1.32 give the impression that the eigenvalue branches which 
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e 1.32: Evolution of 

= t + 0.2(sm7t + 0.5cos27t),7 = 1.1,72 > 7 > 73 



im); 



driver: 



come down cross the eigenvalues, that remain constant in time. This, of course, 
is not possible as the symmetric, tridiagonal operator L cannot have double 
eigenvalues. Instead, a close look demonstrates, that a billard ball collision is 



taking place as shown in Figure |1.33| . It is possible to analyze the interaction 
of Afc and Xk+i in detail by using equation (2^) of Chapter |^ for j = k and 
j = k + 1, together with the asymptotic assumption that |Afc — Xk+i\ is much 
smaller than the distance between any two other eigenvalues, but we do not 
present the details. 




(1.34) 



Figure 1.33: Close look at a "collision" of two eigenvalues 

For A < inf cJess(-^(0)), an interesting quantity to compute is 

{eigenvalues of L{t) that are < A} 



J(A) = lim 

f— >oo 
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Clearly J(A) represents the asymptotic flux of eigenvalues of L{t) across the value 
A. In Chapter ^ we will extend the definition of J to all values of A. 

It is observed numerically that J(A) indeed exists and for 72 > 7 > 73, say, we 



find that J(A) looks as displayed in Figure 1.35 Thus J(A) is constant in the gaps 




Figure 1.35: Numerically observed J(A) in the case 72 > 7 > 73,7 = 1.1 

and indeed we observe more generally that 

(1.36) J(A) = ^ for A in the j*^ gap. 

zvr 

This is a very intriguing fact, reminiscent of the Johnson- Moser gap labelling the- 



orem |JM| in the spectral theory of one-dimensional Schrodinger operators with 
almost periodic potentials (see also the analogous gap labelling theorem for Jacobi 
matrices 0, g). 

Finally we are at the stage where we can describe our analytical results, whose 
goal is to explain the above numerical experiments. In the Toda case with constant 
driving, it was possible, using the exact formulae of inverse scattering, to show that 
the solution of the initial value problem converges as t — > 00 to the binary motion 
if a > 1, and to a quiescent lattice if a < 1. In the present case, where we no longer 
have these formulae, our goals are more modest and we restrict our attention to a 
description of the observed attractor. Our results are the following: 
I. Strongly nonlinear case. 



Here we consider (1.15) in the case of the Toda lattice without any smallness 



restriction on the size of the oscillatory component h of the driver xq. The main 
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result is Theorem |2.38| below in which we show how to compute the normalized 
density of state J(A) through the solution of a linear integral equation, once the 
number and endpoints of the bands in cj(L(oo)) are known. This linear equation, in 
turn, can be solved explicitly via an associated Riemann-Hilbert problem. 

At this stage it is not clear how to relate the number and endpoints of the 
bands to the parameters of the problem a, 7, h. To test Theorem p. 38 in any given 



situation, one reads off the discrete information given by the number and endpoints 
of the bands from the numerical experiment, and then compares the solution of 
the integral equation with the normalized density of states J(A) obtained directly 
from definition (1.33) using the numerically computed eigenvalues of L{t) at large 
times. The numerical and analytical solutions for J(A) agree to very high order: see 
Appendix (Figures p.l2 , C.13 ) for further details. 



The proof of this result proceeds by deriving an equation of motion (see ( p.5[ )) 
for the eigenvalues of a truncated version of L{t) of size N ^ t as t ^ 00. The 
continuum limit of the time average of these equations, leads to the linear integral 



equation (|2.28| ) for J(A). 



II. Weakly nonlinear case. 

Here we consider general F, but the periodic component h is now required to 
be suitably small. From the numerical experiments we see that if /i = 0(e), then as 
t — > 00, Xn{t) converges to an asymptotic state which is a ^-time periodic solution 
Xasymp, n{t) with a^asymp, n (0 = cn + 0(e) for some lattice spacing c. The goal here is 
prove that such time periodic asymptotic states a;asymp,n(*) indeed exist for e small. 
We proceed by linearizing around the particular solution a;asymp,n(i) = cn, n > 0, 
of the equations Xn = F{xn-i — x„) — F(x„ — x„+i), n > 1, and use various tools 
from implicit function theory. 



Our first result (Theorem 3.38| ) is a nonlinear version of the classical linear 



method of separation of variables. For example, in solving the heat equation ut = 
Uxx on a half-line x > with boundary conditions at x = 0, one proceeds by 
expressing the solution combination 



a+{z) e-(*^'+^^^)dz + J a_(z) e^^*^'-^^^) 



of elementary solutions e ('^^i*^^) of the heat equation on the full line, and then 
choosing the parameters a+,a_ to satisfy the boundary condition at x = 0. In the 



nonlinear case (Theorem 3.38|) we show that provided a sufficiently large parameter 
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family of travelling wave solutions of the doubly infinite lattice 



(1.37) 



F{Xn-l - Xn) - F{Xn - Xn+l), -OO < n < OO, 



exist, then (modulo technicalities, see Chapter ^) the parameters can always be 
chosen to produce the desired asymptotic states Xasymp, n(i) of the driven semi- 
infinite problem. 

Thus the problem of the existence of the observed asymptotic states, reduces 
to the problem of constructing sufficiently large parameter families of travelling 
waves of the full lattice equation ( |1.37 ). As we will see in Chapter ^, for 7^ > 
7 > 7fc+i k > 1, we will need 2A;-parameter families of fc-phase travelling waves 
of type ( |1.22D on the full lattice in order to construct the solution of the driven 
lattice observed as t ^ 00 in the numerical experiments. If 7 > 71 (see Section 3.4 ) 
the requirement of travelling wave solutions of ( |1.37 ) trivializes, and Theorem 3.38 
guarantees the existence of the desired asymptotic states 2;asymp,n(i) of the driven 
lattice for sufficiently small h and general real analytic F which are monotone in 
the region of interest, and this explains Figure [Cl6 . 



The next result (Theorem 4.24 ) concerns general F in the case that 71 > 7 > 72. 
Here we show that a 2-parameter family of one-phase travelling wave solutions of 
(1.37) always exist for general F. Together with Theorem 3.35, this implies that 
for 71 > 7 > 72 the desired states Xasymp,n(i) of the driven lattice exist, and this 



explains Figure |C.7| . This 2-parameter family is constructed by deriving an equation 
for the Fourier coefficients of the travelling wave solution, which can be solved by a 
Lyapunov-Schmidt decomposition. The infinite dimensional part does not pose any 
problems (see Lemma 4.19) and the degenerate finite dimensional equations can be 
solved by using certain symmetries of the equation (see Lemma |4.20| ) . 

If we try a similar construction for mo-phase waves, iuq > 1, then we encounter 
in the infinite dimensional part of the Lyapunov-Schmidt decomposition, a small 
divisor problem related to the small divisor problem occurring in | C"VV[| , where pe- 
riodic solutions of the nonlinear wave equation are constructed, and which we hope 
to solve in the near future. In the special case of Toda, however, the family of 
travelling waves can be constructed explicitly. Indeed in our third, and final, result 
(Theorem 5.13| ) we use the integrability of the doubly infinite Toda lattice and show 
how the well-known class of g-gap solutions contains a sufficiently large family of 



travelling waves to apply to Theorem 3.3S and so construct the desired asymptotics 
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states Xasymp,n(*) of the driven lattice for any 7 S R+\{7i, 72, • • •}. 



Finally we want to pose four open problems, which are connected to our inves- 
tigations, some of which were mentioned above. 

(1) The "critical shock" phenomena. 

As discussed in the very beginning of the introduction, Holian and Straub 
have numerically discovered a critical shock strength acrit(^) in the case of 
constant driving velocity xoit) = 2at. For a < acrit(-^) the lattice comes to 
rest behind the shock front as t ^ 00, whereas for a > acritl-P") the particles 
of the lattice will execute binary oscillations as t ^ 00. So far this result 
has been analytically explained in the case of the Toda lattice {F(x) = e^) 
(cf [ pFM |, [VDO|) and can be seen to be absent for linear force functions by 



explicit calculation. 

The question is to find general conditions on the force F for which one can 
prove the existence of a critical shock strength. 

(2) Existence of multi-phase travelling waves. 

Let F, c satisfy the general assumptions (cf Section 3.2.2| ) and let 7 G M satisfy 
ly/F'{-c) > 7 > j^^/F'{-c) for some kGN. 

Does there exist a smooth 2A;-real parameter family of solutions 

(1.38) Xniq)it) = cn + Xk{q){n(3i{q) +jt,.. . ,n(3k{q) + kjt), 
for q S M^'^jq small, of the equation 

(1.39) Xnit) = F{xn~i - x„) - F{xn - x„+i), n G Z, 

where Xk{q) is for each q a function periodic in its k arguments, XA:(0) = and 
DqXkiO) has maximal rank 2k'! Note that these solutions exist in the case of 



the Toda lattice and are given by ( 5.10 ). For general force functions the work 



of Craig and Wayne (| CW|| ) indicates that it might be necessary to aim for a 



slightly weaker result, namely, that the smooth family of functions Xniq){t) of 
the form (1.38) are solutions of ( [L.39 ) only for a Cantor set in the parameter 



space q G R'^'^, which has almost full measure. 
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(3) The case a > Ccrit- 

In the paper we always assume that the driver is of the form XQ^t) = 2at+h{"ft) 
with a < Ccrit- In the case of the Toda lattice we have also conducted some 
experiments for a > acritiToda) = 1. We have made the following observation: 
for h small, the limiting operator L{t),t — > oo seems to have infinitly many 
gaps and again we obtain a version of gap labelling. In fact, for all the gaps 
we have observed that one can write the numerically determined integrated 
spectral density J(A) (cf |1.34| ) in the form 

(1.40) J(A) = — , j, k £ Z, A lies in a gap, 

zvr 

where 7 denotes the frequency of the driver and uj is given by the frequency 
of the time-asymptotic oscillations, which are observed in the case that the 
driver has constant speed xo{t) = 2at, a > 1. 

Corresponding to our results in Chapter |5|, we ask whether it is possible to 
construct solutions of the driven semi-infinite Toda lattice with driver xo{t) = 
2at + h{-^t), such that the spectrum of the corresponding Lax-operator has 
infinitely many gaps and bands. 

(4) Connection to the initial value problem. 

All of our results were motivated by the numerically observed long-time be- 
havior of a certain initial value problem (with shock initial data). However, 
so far we are not able to prove from first principles, that the solution of the 
initial value problem actually converges as t ^ cxd to one of the asymptotic 
states described in Chapters |2[|5|. This basic problem remains open and, alas, 
seems far from a resolution. 
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Chapter 2 



An asymptotic calculation in 
the strongly nonlinear case 

2.1 The evolution equations 

We recall from the Introduction the Flaschka variables (see (|1.12| )) 

(2.1) «- = -Y' = ^e(— n = 0,1,2,..., 

On — > a, bn ^ ^ as n — > +00, 

and we note that the function ao(^) is the given time-periodic forcing function. For 
the semi-infinite Toda chain with F{x) = e^, equation ( pTl] ) reduces to the perturbed 
Lax pair equation 

(2.2) ^ + LB-BL = -p{t)P, 
where L is the tridiagonal operator (cf (|0^)) 

^ ai 61 ^ 

hi 02 62 

L = 

62 as 63 

V ■■• ■•• ■■• ) 
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B is the antisymmetric tridiagonal operator given by 



B = 



/ 61 
-61 62 
-62 

V 

P is the rank-one matrix given by: 



63 







if i = j = 1 
otherwise 



and p{t) is the function: 



(2.3) 



p{t) = 2bl{t) = 2blit)-ai{t), 



A. 
dt 



The matrices L and B are semi-infinite. We truncate the chain at some particle 
of very large index N, and work with the truncated finite matrices Ljv and B^. 
The disturbance in the chain caused by the truncation, travels essentially with 
finite velocity. Only exponentially small effects display infinite speed. The bulk of 
the chain essentially does not feel the truncation until a time t = 0{N). Thus, we 
expect that the finite system is a good approximation to the full semi-infinite system 
in the space-time region 1 « t « N and n « N. In what follows, when we take 
the limit as t —> 00, we always understand t ^ og, N co, ^^0. 

We remark that on aessiL{t)) = aessiL{0)), by standard spectral methods, the 
matrix L^lt) has a set of discrete eigenvalues tightly packed at densities of order 
0{N). On the other hand, we expect that the discrete spectrum of L{t), which 
emerges from cjess(^(0)) as described in the introduction, is well approximated by 
the eigenvalues of Li\f{t) which lie outside cJess(-^(0)). This is because the associated 
eigenvectors are typically exponentially decreasing in n and hence do not feel the 
truncation at N. 

A word of explanation: Typically an eigenvalue of Li\f{t) starts off as an eigen- 
value of L7v(0) lying in aess{L{0)). As t increases, the eigenvalue moves with velocity 
0{j^) until it emerges from aessiL{0)). It is only after this point that the motion of 
the eigenvalue becomes relevant to the evolution of the discrete spectrum of L(t). 

Our strategy is to derive evolution equations for 



(a) the eigenvalues of the truncated matrix L^, 
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(b) the first entry fj of the j^^ eigenvector of L]\f (j = 1, . . . , A'') when it is 
normalized to have Euclidean length equal to one. 

It is well known that the set {Xj, determines the tridiagonal matrix Ljv- 

Theorem 2.4 The evolution of the Xj 's and fj 's is given by: 

\ ^|ln(-A,) = A,-ao(t)+Ef!i^, J = l,...,iV, 

(2.5) <^ 

\ Jj pi 

where p = 26g(f) = — X^i^i ^i- initial values Aj(0) are the eigenvalues of L]\f at 
t = while the initial values Aj(0) are given by 

(2.6) A,(0) = -262(0)/2(0). 

Proof : Let A be the diagonal matrix of the eigenvalues Xj of L^r and let ^ 
be the orthogonal matrix whose j^^ column is the normalized eigenvector of Ljv 
corresponding to the eigenvalue Xj. We have: 

(2.7) Ln"^ = ^-A, 
and we define the matrix ^ by: 

(2.8) ^>=^'-SAr^'. 



Utilizing equations (2/7) and (|2.8D and (|2.2| ), we easily calculate: 

(2.9) Ljv^> - ^>A = ^A + pP^-. 
We now define the matrix A = {aij) by 

(2.10) ^ = ^(-1$ = 
We calculate 



- Bn^) + [^^ - Bj^)^ 
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Thus 

(2.11) A + A'^ = 0; 
i.e. A is antisymmetric. Using ( p.9| ) we obtain 

(2.12) Ljv^>-^>A = Ljv^'yl-^'^A = ^'(AA-^A). 
Comparing ( |2.9D with ( p. 12 ) we obtain easily: 

(2.13) A= [A,A] -p^'^P*. 

Let = (/i, /2, . . . , /at) be the first row of Then ^''^P^' = //^. We insert this 
relation in ( |2.13| ), 

(2.14) A=[A,A]-p//^. 
Equating the diagonal elements on both sides we obtain: 

N 



(2.15) 



A,- 



This proves the second relation in Theorem 2.4. Furthermore we note that the first 
components fj of the eigenvectors of the tridiagonal matrix do not vanish and 
we conclude by ( 2.15| ) that —Xj > 0. Hence the first relation in Theorem 2.4 is well 
defined. 

Off the diagonal in ( p.l4| ) we have AjOij — aijXj = pfifj. On the other hand 
an = by pTTl] ). Thus 

_ pfifj 



(2.16) 



0. 



when i 7^ j, 



We now calculate the evolution of fj . By (^) : 

^ = $ + Bn"^ = ^A + Bn"^- 

Specializing to the first row we obtain f'^ = f'^A + Pjj^^, where Bji^ is the first 
row of Pat. This implies 

= fA+i^A)n,-aif = fA + fA-fau 



or, taking transposes 
(2.17) 



f = {A-aiI-A)f. 



23 



But, by (|2ll ), 
(2.18) 



A = pFDF, 



where F is the invertible diagonal matrix with entries fi, ■ ■ ■ , /n and D is the matrix 
^ ' with zero entries on the diagonal. Thus / = {A — ail — pFDF)f , and hence 



(2.19) 

Now note that 
(i) pFf = p 



(ii) F-i/ 



fl 

V '■ J 

1 

V '■■ J 



F-'f = F-'Af - aiF"V - pDFf. 



A2 

V J 



F^^A = AF- 



by (EH), 



Substituting in (|1|) we obtain 
(2.20) 



The evolution of the Ai's in (2^) is finally obtained by eliminating fj between ( p. 20 ) 
and ( |2.15 ) and using the expression for p given in (^.3|). 



O 



2.2 The continuum limit of eigenvalue dynamics 

The results of the numerical experiments described in the Introduction, (cf Figures 



1.2£ -( 1.32 lead us to consider the fiux of eigenvalues of the matrix Ljv across a value 
A. Noting that the eigenvalues of L^r move toward lower values (Aj = —pff < 0), 
we define the eigenvalue fiux at A averaged over a time interval {t,t + T) by 

(2.21) JtT{X) = ^ card {i : \i{t + T) < \ < Xi{t) , i = l,2...N}. 

We pose the following ansatz. 
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Ansatz 2.22 There exists a continuous, almost everywhere continuously differen- 
tiahle function J(A) such that 



(2.23) Jt,T(A) ^ J(A), when T,t,N^(x> subject to T<t<^N. 
When A < inf aess{L), it is clearly true that 

(2 24) J(A) r eigenvalues of L{t) that are smaller than A 

as defined in the Introduction. The net gain in eigenvalues of of an interval (A, A) 
over a long time T is given asymptotically by [J(A)— J(A)]T. Dividing by (A— A)T and 
letting A ^ A we obtain that the asymptotic rate of increase in eigenvalue concen- 
tration at A is given by ^ = J'(A); thus, the difference in eigenvalue concentration 
at A between times t and zero is asymptotically tJ'{X). When A < inf aess{L), 
necessarily J' (A) > since there is no original eigenvalue concentration at A. On 
the other hand J'(A) can take negative values when A G aess{L). 

We will now use the function J(A) and some assumptions based on numerical 
observations to derive the continuum limit of the eigenvalue evolution equations 
(|2.5|). We begin by averaging the system (2^) of equations (j = 1, . . . N) over the 
time interval {t,t + T) to obtain 

1 XAt + T) 
(2.25) T^ln - 



2T A,(t) 



i-t+T 1 rt+T ^ JL rt+T (iXiif) 



-J^ X^{t')dt' - - ao{t')dt' + -J2j^ 



W)-Xi{t') 



l,...iV. 



»7^i 

Let A satisfy J'{X) ^ 0, and let j = j{t) be such that in the asymptotic limit 
I <^ T <^ t <^ N (note that we require T <^ t, not just T < t as in ( p.23| )) the 
following is true, 

(2.26) ^j{t){t') X for any t' that satisfies t<t'<t + T. 

In practical terms, this means that we expect eigenvalues Xj to stay close to the 
value A throughout the time interval [t, t + T]. This fact is clearly borne out in the 
results of numerical experiments as long as J'(A) ^ 0. 

We make two more simplifying assumptions when J'(A) ^ that are again 
justified by numerical experiments: 
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(a) The left hand side of (p.25 ) is negligible. 

(b) The "singular contribution" in the sum of the right hand side correspond- 
ing to indices i that are close to j is also negligible. In practical terms 
we interpret this to mean that the limiting integral kernel is the Hilbert 
transform. 



Under these conditions we can take the limit 1 <C T < t < iV in (|]2|)-(|]2|). 
Theorem 2.27 (Continuum Limit of ^.2^ )-( ^.2(\ )). Under Ansatz ^^^^ and un- 



der the further assumption described above we have: 

no 

Jifi) 



(2.28) 



J'(A) / O^A-<oo> - / Y^'^^ = 0' 

J A — fi 



where < aQ > is the mean value of the periodic driver ao{t), and as usual the double 
bars on the integral indicate that the principal value is taken. 



Proof : By the assumptions, and by ( 2.26| ) the only thing to be shown is that 



the sum in ( 2.25 ) tends to the integral in ( 2.2^ ). If we partition the eigenvalue axis 



into a set of infinitesimal intervals and if (fi — dfi, fx) is such an interval then the 
contribution — ( - — —] should arise in as many terms of the sum in (|2.26) as there 



are eigenvalues that cross the value /i during the time interval (t, t+T). This number 

oo 

is asymptotically TJ{ij,). The sum in ( |2.26| ) therefore tends to — ^ J{p)dfi ^ 

J A — /i 



2.3 The asymptotic spectral density of L 



N 



We now consider the problem of determining J(A). Our solution is partial in the 
sense that we can calculate the function J(A) if we are given the set {A : J' (A) ^ 0}. 
Numerical calculations (see Figures [L31| , |L3|, [OsI) show that this set is a finite 



union of intervals. Thus, we are assuming knowledge of a finite set of numbers 
that are in principle determined by the fluctuating part of the periodic driver ao(t). 
Determining these numbers is, unfortunately, the part of the problem that we have 
not yet been able to solve. 

We proceed to give some basic definitions. 



26 



Definition 2.29 Let the points po < qi < pi < q2 < P2 < ■ ■ ■ < Pg < Qg+i be 
given. These are 2g + 2 points in all. We define the set o/ bands B, where we have 
J'{X) / 0, 

(2.30) B =[po,qi]U[pi,q2]U---U\pg,qg+i], 
and the set o/gaps G, where J'(A) = 0, by 

(2.31) Gk = {qk,Pk), k = 1,2,... g; G = Ul^^Gk 
(cf Figure 2.3(\ below). We then define the hyperelliptic curve 

(2.32) RiX) = {n^^o(A - pk){X - gfc+i)}^/' 

with branch cuts along the set B and sign determination such that R{\) > when 
A — > + oo. Finally, we define the polynomial 



(2.33) 

where the ai's, i = 1, 
(2.34) 

(2.35) 



9+1 



P(A) = 5](A-aO ,(TiGR, 



i=l 



. g + \, are uniquely determined by the g + 1 relations. 
P{\) 



RiX) 



dX = ,k = l,2,...g. 



B 



R{X) 



dX = 0. 




Po qi Pi 



Figure 2.36: The band - gap structure of the spectrum 
We observe that the integrals in ( 2.34| ) and ( |2.35| ) can be easily understood 



as contour integrals on the Riemann surface associated with R. The contour in 



P(A) 



( |2.35D can be replaced by a circle of (large) radius. We then obtain -^^j-^^ 



1 + 



(2.37) 



as A — > cxD, which implies through an easy asymptotic calculation that 

9+1 9 
2^1'^* = '^iPi + Qi+l)- 
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Theorem 2.38 Let J(A) be a continuous function supported on the set B L) G, 

differ entiahle at all points, except possibly the boundary points of B, satisfying 



(2.39) 

(2.40) 
(2.41) 



A- < ao > - ^ 4^dfi = X£B, 
J A - /i 

— oo 

J(A) = Cfc = const., X £ Gk, k = 1 . . . , g, 
J(A) = \^ BUG, 



where < uq > is a constant. Then J{X) + iHJ{X) is the limiting value as z ^ X + iO, 
of the analytic function 



(2.42) 

Precisely: 
(2.43) 

(2.44) 



oo 

J(A) = ilm/(l 



dfi, 



The endpoints of B satisfy the compatibility condition 

R{X)' 



(2.45) J (^1 - ll^dA ) =< ao > 



99 + 1 



Remark 2.46 on condition ( 2.4(\ ). 

In the above Theorem 2. 3S^ we have not specified the value Ck which the function J(A) 
obtains in the k-th gap. In fact, Ck will be determined by all the other conditions 



of Theorem 2.3i . However, as remarked in the Introduction ( 1.36 ), one observes 



in numerical experiments that c^ should equal kj/2TT. Figures C.li and CIS in 
Appendix demonstrate that the solution J(A) of the integral equation ( 2.3!^ - 



(2.41) indeed satisfies this additional relation. 



Proof : The differentiability properties of J(A) in the interior of B as well as its 
constancy on the G^s is immediately obvious as soon as one sees that R is pure 
imaginary in the interior of B, and pure real elsewhere. The function J(A) is clearly 
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continuous at the endpoints of each Gk and at qg+i- Also J{po) = by (|2.35| ) and 
consequently J(A) = when A < po- Condition ( 2.39|) follows from (|2.34| ), ( ^.44 ) 



and ( p.45 ) in a straightforward way, using once again the pure real/pure imaginary 



structure of R. 

The function J(A) constructed is unique. Indeed, if by 6{X) we denote the 
difference of two solutions of ( 2.3g| )-( |234T ), then d{X) satisfies the equations 



= 0, when X £ B, 

(2.47) ^ 6{X) = c'f^ = const., when A G Gk, n = 1, 

= 0, when X^BUG, 




and its derivative 6'{X) satisfies 
(2.48) 




0, when Xe B, 
0, when X ^ B. 



By ( |2.48| ) and the third relation in ( 2.47] ) the derivative 6'{X) of 6{X) is identically 



zero; since 6{X) is compactly supported we also have (5(A) = 0. 

O 

The above derivation of the equation for J(A) has been obtained under a variety 
of assumptions. A fully rigorous derivation still eludes us. At this stage however 
the justification lies in the comparison with experiments as described e.g. in Figures 
aT2| and of Appendix g 
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Chapter 3 



A boundary matching 
technique 



3.1 Introduction 



In this chapter we construct - time periodic solutions of 

(3.1) Xn{t) = F{Xn-lit)-Xn{t))-F{Xn{t)-Xn+l{t)), n > 1, 

with 

(3.2) xo(t) := e J] bme'^^\ Yl l^-l = 1' 

mGZ JTiGZ 

satisfying 

(3.3) Xn{t) = cn + 0{e). 

Throughout this chapter we assume that F is real analytic and monoton increasing 
on an open interval. We will construct solutions satisfying (|3.1|)-( |3^ ) for any c € M, 
such that — c lies in the interval. Not all values of c, however, can be observed as 
the spacing of an asymptotic state of the driven lattice, described by the initial 



boundary value problem I.IJ. To see this we look at the Toda lattice. For a < 1, 
in the case where h = 0, the solution Xn{t) converges as t ^ oo to Xasymp,n(0 = cn, 
with spacing c = —2 ln(l + a). Thus the values of the spacing c that can be observed 
by driving the Toda lattice with constant velocity 2o, lie between — In 4 and 0. 
For a > 1, as we know, the solution of the driven lattice does not converge to a 
quiescent state Xasymp,n(i) = cn and, in particular, the values c = — 21n(l + a), for 
a > 1 cannot be observed in this experiment. 
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Expand Xnit) in a Fourier series, 

(3.4) Xnit) =cn+J2 rn)e'"''^\ for n > 0. 

mGZ 

For n = we have a{0,m) = ebm- Expanding F in a power series at — c, we obtain 

= F'{—c) ^ [a{n — 1, m) — 2a(n, m) + a(n + 1, m)]e*'^™'* + higher order terms . 
Equation (^) with (|3.2| ) is equivalent to 

(3.5) (Lma(-, 'm.))(?i-) + VF(a)(n, m) + a(0, m)5i^„ = 0, for all m £ Z, re > 1, 
where for m G Z the linear operators acting on the re-variable are given by 

/ 5n,, 1 \ 

(?re7)^ 



(3.6) 



1 5« 



1 



1 



V 



-2 + 



F'i-cY 



J 



W{a) contains all terms of higher order and (5i „ denotes the Kronecker symbol. 

We note that for e = equation ( |3.5| ) is solved by a = 0. However, we cannot 
apply the implicit function theorem to obtain solutions of equation ( p.5| ) for e 7^ 0, 
because the linearized operator, 



(3.7) 



is not invertible. Indeed, the spectrum of the operator Lm acting on I2 sequences 
is given by a{Lm) = [6m — 2,(5^ + 2]. This implies that G cr{Lm) for all m G Z 
satisfying < piJ^J^^ < 4. Denote 



(3.8) 



mo := maxirei, E Z : < X , n < 4}, 

F'(-c) 



then the multiplicity of in the £2 spectrum of L is 27reo + 1. However, due to the 
simple form of the operator L^n we are able to " invert " the operator explicitly. In 
fact, consider the important case < pi"^^-^ < 4. For a given vector {yn)n>i and 
given ui, the vector (?i„)„>2 solves the equation 

LmU = y, 
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if and only if 



1 



sin P„ 
1 

sin (3„ 



n-l 



ui sinnPrn + ^yk sin(n - k)(3r, 
k=l 

sin [ui + ^yk cos k/Sm - cos n/?m ^ sin A:/3„ 



n-l 



fc=l 



fc=l 



(3.9) 



/3« 



sgn (m) arccos 



Note, that the value of Un is independent of the sign of (3m- We will justify the 
particular choice we have made in Section |3]^ (see ( p.35| )) below. 

The following observation will prove to be useful. Suppose (yn)n>i decays 
exponentially, i.e. there exists a cr > 0, such that sup„>i \yne'^^\ < oo, then 
sup^>]^ I'Une'^"! < oo, provided the following two relations hold. 



(3.10) 

and 
(3.11) 



ui + '^Vk cos k(3m = 



k=l 



Vk sin kPm = 0. 



k=l 



Equation ( p^ ) can always be satisfied by an appropriate choice of ui, whereas 
equation ( |3.11| ) is a condition on the sequence (yn)n>i- Therefore the operator 
acts 1-1 on spaces of exponentially decaying sequences and the range has codimen- 
sion 1. Furthermore a simple calculation shows that the inverse operator acts on 
the range as a bounded operator with respect to the corresponding exponentially 
weighted supremum norms. Still we cannot apply a standard implicit function the- 
orem to obtain solutions of equation (|3.5|). Nevertheless, proceeding formally, we 



transform equation (3.5) into a fixed-point equation. 
(3.12) a{n,m) = -L:;^^[W{a){n,m) + ebmS 



l,nj 



As W is of higher order we can in principle apply a Banach fixed-point argument 
to obtain a solution of ( 3.12|) as long as is a bounded operator. We have seen 



above that this can be achieved, if condition (3.11) is satisfied, i.e. 



(3.13) 



ebm sin + ^ W{a){k, m) sin kf3m = 0, for < 
k=i 



(7m)^ 
F'i-c) 



< 4. 
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Equation ( ^.13 ) indicates that we will be able to solve equation ( p. 12 ) for sufficiently 



small e in a space of sequences decaying exponentially in n, only if the Fourier 
coefficients of the driver bm take on a special value for those m £ Z satisfying 

This observation is consistent with the linear case where we conclude from for- 
mulae (|1.23| ) - ( 1.25 ) that the solutions decay exponentially in n, only if 6^ = for 



< 1^ < 4. The hnear case also suggests that we should add multiphase waves 
in order to obtain solutions of equations ( |3.1D ,( |3^ for general driving functions. 
This leads to the following ansatz for a{n,m): 

(3.14) a{n, m) = u{n, m) + v{n, m) + {ebm — u{0, m) — v{0, m)), n > 0, m G Z, 
where u denotes the travelling wave part and v corresponds to the exponentially 



decaying modes. Note that ( 3.14 ) implies a(0, TTi) — ^byyi-i for all Tfi G Z. 



Definition 3.15 We will refer to the Fourier modes m with < pi^^-^ ^ 4, or 
equivalently \m\ < rriQ, as resonant Fourier modes. On the other hand we say 
that a frequency 7 G is resonant if pi^^^ = 4 for some m G Z. 

The present chapter is organized as follows. We begin Section |3.2| by deriving 



equations for the sequences of Fourier coefficients u{n, m) and v{n, m) (given by (3.4) 
and ( p. 14 ) above), which are sufficient to prove that the corresponding functions 
Xn{t) solve (3.1), ( 13.21 ). These equations, which are given in Lemma p. 31 below, can 



be made rather explicit because of the assumption, that the force function F can 
locally be expanded in a power series and therefore we will obtain good estimates on 
the higher order terms by carefully choosing the norm on the sequences of Fourier 
coefficients. In the notation of Lemma p. 31 , these equations can be described as 
follows. 

(1) is an equation for n, which is satisfied by the Fourier coefficients of solutions 
Xn{t;u) := X^mGZ ""(^^ the doubly infinite lattice. 

(2) is an equation for v, depending on u, which guarantees that u + v corresponds 
to a solution of the semi-infinite lattice. 

(3) represents the boundary condition by requiring ebm — u{0,m) — v{0,m) = 0, 
for m 7^ 0. The case m = is special; we do not have to require ebo — 
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m(0, 0) — v{0, 0) = for the reason that solutions of ( |3.1| ) are invariant under 
translations x„ — > x„ + const. 



We then proceed in Section 3.3 to prove the basic result (Theorem |3.38| ) of this 



chapter. Assume 7 is non resonant (see Definition 3.15| , then for (small) e and 



for given (small) travelling wave solutions u of the doubly infinite lattice, we can 
construct sequences v satisfying equation (2) above and solving equation (3) for all 
non resonant Fourier modes m (compare with Definition 3.15| ), i.e. for those m G Z 



satisfying \m\ > ttiq. Furthermore suppose that u is given as a function of a 
parameter q. Then we will show that the resulting v is a function of q and e. 
Note that this statement is needed in order to ensure that the remaining equations 
of (3) (for resonant Fourier modes) can be solved by constructing a sufficiently large 
parameter family of travelling wave solutions u{q), and then applying a standard 
implicit function theorem. 



The proof of this basic result (Theorem 3.38| ) rests on a Banach fixed-point 
argument. The equation for v takes the form 

(3.16) {Lv){n, m) + W{u, v){n, m) + t>(0, m)5i^n = 0, for n > 1, m G Z. 



L = @rn&^m was defined in (13. q) and (^) and W{u,v) denotes the higher order 
terms. We turn ( |3.16| ) formally into a fixed-point equation, 

(3.17) v{-,m) = -L:^ {W{u, v){-,m) + t;(0, m)(5i,.) . 

Denote 

Sa ■■= { {yn)n>l ■ SUp IVnG^'^l <Oo \ . 
[ n>l J 

As indicated above we will be able to prove the following results on the invertibility 



of Lm by explicit calculation (see proof of Theorem 3.38) 



For < |m| < niQ and o" > 0, the linear operator maps onto {y £ ■ 
J2k>i Vk sin(A;/3m) = 0}. (The quantities j3m were defined in (p.Op. The inverse 
operator acting on the range is bounded with respect to the corresponding 
norms. 

There exist weights o" > 0, such that the operators Lm : ^ are bijective 
and have a bounded inverse for all \m\ > mo. 
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• Again the case m = is somewhat special as the Green's function of the 
operator Lq grows linearly and we will have to use the special structure of 



W{u, v) in order to define a bounded inverse. See the proof of Theorem 3.38 
below for more details. 

Although n(n, m) does not decay in n, we will nevertheless see by explicit calcu- 
lation that W{u, v) decays exponentially in n. This makes it possible to prove the 



exitence of a solution of ( 3.17 ) by a Banach fixed-point argument. For \m\ > ttiq 
we can choose v{0,m) = ebm — u{0,m) and hence satisfy the boundary condition 
as described in equation (3) above, whereas in the case < |m| < niQ the choice 
of w(0, m) is determined by the condition that W{u,v){-,m) + v{0,m)6i,. has to lie 
in the range of Lm- A small technical problem arises when proving the smooth de- 
pendence of V on the parameters. It will turn out that the travelling wave solutions 
constructed in the subsequent chapters depend smoothly on q, but §^(ra, m) grows 



linearly in n. Therefore does not lie in a space which is suitable for our calcula- 



dq 

tions. We will verify the smooth dependence of f on g and e explicitly by applying 
a Banach fixed-point argument to the partial derivatives in the appropriate spaces. 



In Section 3.4 we show that the results of Sections and suffice to construct 
periodic solutions of the driven lattice in the case that mo = 0, which corresponds 
to high driving frequencies. 



3.2 The equation for the Fourier coefficients 

In this section we introduce norms, which are suitable for the sequences of Fourier 
coefficients, and prove some of their basic properties. Then the general assump- 
tions on the force function F and on the driver will be stated precisely. Using the 
assumptions on F we derive estimates on the nonlinear terms which allow us to 
give conditions on the Fourier coefficients which are sufficient for proving that the 
correponding functions x„(i) given by (|3.4D and ( |3.14|) solve equations ( ^T| ) and 

dH). 

3.2.1 Sequence spaces 

The choice for the norms on the sequences of Fourier coefficients u{n, m) and v{n, m) 
(compare with equation ( p. 14] )) is motivated by the following observations. The 



35 



nonlinear terms of the force function make it necessary to take convolutions with 



respect to the m- variable (see Section 3.2.3| below). Therefore we choose an ^i-norm 
for m. In fact we use a weighted £i-norm in order to control the regularity of the 
solution. Furthermore the weight function has to satisfy some additional conditions 
to insure that the norm is still compatible with respect to convolution (see Definition 
3.1^ of admissible weight functions) . For the n- variable a supremum norm with an 
exponential weight is chosen which is suitable for inverting the linearized operators 



Definition 3.18 A map w : Z 



is said to be an admissible weight function, if 



(3.19) 

and 
(3.20) 



w{m) > 1, for all m G Z, 



w{m) < w{m — n)w{n), for all m,n £ Z. 



Definition 3.21 Let w be an admissible weight function. We denote 



-l,w 



u:Z^C 



w{m)\u{m)\ < oo > , 



mGZ 



with the corresponding norm 



w{m)\u{m)\. 



mG2 



Note that is a Banach space which lies in ii by condition ( 3.1S| ). The in- 
equality (|]2|) insures that the ii^ui-noim is submultiplicative with respect to con- 
volution (see Proposition p. 23 below). It also implies that the weight function 
can not grow faster than exponentially. Indeed, it is easy to prove that w{m) < 
w{0) (max(t(;(l), ti;(— 1))) . We shall be interested in three types of weight func- 



tions which will all satisfy the conditions specified in Definition 3.18. 

(i) Vm G Z : w{m) := 1 . 

(ii) \/m£Z: w{m) := (1 + |m|)^, for /? > . 

(iii) Vm G Z : w{m) := e^l^l, for f3>0. 

Finally note that the product of two admissible weight functions is again an admis- 
sible weight function. 
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Definition 3.22 Let w be an admissible weight function and let o" G R. Then 
^a,w := < u : No X Z ^ 



supe™|u(n,-)| G^i 
> ?i>0 



mGZ ) 



with the corresponding norm 

w{m) sup 6°^" I n(n,m) I . 



I c^^w 

n>0 



It is easy to check that Ca,w are Banach spaces and that C^i^w C £0-2, w foi^ '^i ^ cr2- 
In the following proposition we recall some simple properties of the convolution of 
sequences, which is defined by {u * v){m) := J2iez'^i^ ~ 0^(0- Furthermore we 
provide estimates on the convolution in terms of the norms defined above. 

Proposition 3.23 Let c £ C,ui,U2,U3,U4 £ ii, then 

(i) ui * U2 £ l\ and \\ui * U2\\e-^ < ||ni||£j|n2||£^. 

(ii) Ui * U2 = U2 * Ui. 

(Hi) ui * {u2 * U3) = (ui * U2) * U3. 

(iv) ui * {u2 + CU3) = ui* U2 + c{ui * U3). 



(v) If ui and U2 satisfy the reality condition (i.e. Vm G Z : Ui{m) = Ui{—m);i = 
1,2), then ui * U2 also satisfies the reality condition. 

(vi) If for all m £ Z : \ui{m)\ < U2{m) and \u3{m)\ <U4^{m), then \{ui*U3){m)\ < 
{u2 * Uii){m) for all m £ Z . 

(vii) Convolution respects the ii-a, norm, i.e. let vi,V2 £ ii^w then vi * V2 £ ii,w 
and \\vi *W2||^i,,„ < ||'yi|ki,„lb2||£i,„ . 

(via) Let (Ti,(T2 £ K and u £ Caj^^w^v £ Ca^.w CL^d define their m - convolution y 
by y{n,m) := Y.i(^zu{n,m - l)v{n,l). Then y £ Cai+a2,w and \\y\\a^+a2,w < 

\\u\\a-i ,w\\v\\a2 ,w ■ 

Proof : Properties (i)-(vi) are standard. In order to show (vii) we note that 
the inequality ( |3.2[lD implies 

^'^w{m)\{ui * U2){m)\ < w{m — n)\ui{m — n)\w{n)\u2{n)\ 

m m,n 

= * (^^1^21)11^1, 

and by (i) this is all we need. Property (viii) is a consequence of (vi) and (vii). 



37 



o 



3.2.2 The general assumptions 



Recall the notation which was introduced in equations ( p.lD and (|3.2| ). We now state 



the assumptions on the force function F, the frequency 7 and the Fourier coefficients 
(ftm)m6z of the driver. 

The general assumptions. 

(1) F : K ^ M is real analytic in a neighborhood of — c, c G R, and 



(3) (6m)mez £ ^i,w for some admissible weight function w and ||6m||£i ^ = ^■ 
Remark: 

• We are looking for a solution of the type = cn+0(e). Therefore F{xn-i — 
Xn) = F{—c + 0(e)). Condition (1) will allow us to expand F{xn-i — x„) — 
F{xn — Xn+i) in a power series where the linear term does not vanish. 

• It will be shown that the exceptional set of resonant frequencies (see Definition 
3. 151 ) {7 : = 4 for some m G Z} consists of precisely those frequencies for 
which the number of phases in the travelling wave solution described above 
changes. In the case of the Toda lattice these are also the frequencies for 
which the number of gaps in the spectrum of the corresponding Lax operator 
at t = cxD changes. 

• It turns out that the weighted spaces li^yj are well suited to proving that the 
regularity of the solution is comparable to the regularity of the driver. 

3.2.3 The nonlinear terms 

The force function F is assumed to be a real analytic function at — c (see general 
assumptions above) and we can define for all A; > 0, 



(3.24) 




(2) 7GKn{7:|^ 



4 for some m G Z}. 



(3.25) 
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By Pf,c we denote the minimum of 1 and the radius of convergence of the power 
series Y^^=o l^i-^ ~^ ■ Recall that ai 7^ by the general assumptions. Therefore 
we obtain the following estimates by standard arguments for power series. 

Proposition 3.26 There exists a constant Cf,c, such that for all y, \y\ < 

1 00 I I 
00 I I 

We now define the higher order terms of the equations for the Fourier coefficients 
as formal power series. The convergence of these series and various differentiability 
properties will be discussed in the subsequent proposition. In order to see that the 
following expressions indeed represent the higher order terms of the equation, one 



may look at Lemma 3.31 below. 

If ti = u{n,m) G Ctj^wi we use u{n, •)*'^ to denote the k -th m-convolution of n, 
that is 

(3.27) u{n,-Y\m)= ^ n(n, /i) • . . . • n(n, Z^). 

iiH |-/fe=m 

Definition 3.28 For a > 0,u,v £ C^^w <ind n>l, denote 
Au{n,m) := u{n — l,m) — u{n,m). 

Wiu)in,m) := — ^ ^ ((A'u)(n, (A'u)(n + 1, (m). 

"1 k=2 

1 00 

Y{u,v){n,m) := _ ^ |: ((A(n + ^;))(n, O*'^ - (An)(n, •)*'^) (m). 
W{u,v){n,m) := Y{u,v){n,m) — Y{u,v){n + l,m). 



For n = and for all m £ Z set Au{0,m) := W{u){0,m) := Y{u,v){0,m) := 
W{u,v){0,m) := 0. 

Proposition 3.29 There exists a constant Cf,c, such that for all < a < l,u £ 

£o,to (ind V e £o-,«> with \\u\\o^w,\\v\\(7,w < the following is true. The series in 
the definition of W and Y converge absolutely with W{u) £ Co^w,Y(u,v) £ Ca,w 
Furthermore 
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(i) \\w{u)h,^<CFAMl.,- 

(iii) The map Fi : {v ^ Ca^w '■ \\v\\(t,w < ^-f^} ^ ^a,w,v Y{u,v) is and the 
derivatives satisfy the following estimates 

Vx e Ca,w : \\{D^Y){u,v)x\\a,w < Cf,c niax(||n||o,^, ||t'||o,u;)||2;||(7,u.- 

yXl,X2£Ca,w ■■ \\iDlY){u,v)[xi,X2]\\a,w < CF,c\\^l\\a,w\\x2\\a,w 

(iv) The map F2 : {u ^ £o,to : ll'wllo.w < ~^ ^a,w,u \—>- Y{u,v) is with 
derivative 

1 00 fc— 1 I I, \ 
(3.30) DuY{u, v)x = —y^yi \{k- /)(A^/)*(^-^-i) * (Av)*^ * Ax. 



DuY{u, v) as given in equation ( \3. 3(\ ) can he regarded as a hounded linear 
operator from Ccji into C^+a' ,w for o"' G M and the corresponding operator 
norm is hounded hy (1 + e*^' )Ci;'^c iiiax(||n||o,«;, ||'v||(j,ti>)- 

(v) Let 0"' > and fix x ^ The map 

F'i '■ ^a,w l^a-cj' ^w^v ^ {DuY){u,v)x is and the derivative satisfies the 
estimate 

,c|| ^ II — (t',UI II -^ll (T,UI • 

Remark: The differentiability properties (iii)-(v) will not be used in the present 
section, but they are needed in Section p.3| when we prove differentiability of the 
solution of the fixed-point equation with respect to certain parameters (compare 
with the proof of Theorem 3.38| ). 



Proof : (i) We begin by remarking that u G £o,to implies An G £o,to and 
|Ati||o,«, < 2||ti||o,«,. By Proposition 3 .231 (viii) and Proposition 3.26 it is easy to 



see, that 



ir(n)||o,. < 2— ^^IIAnllg,, 

l"l| fc^2 
< 2C'F,e(2||u||o,«,)^ 



(ii) In this case one has to evaluate 

iAiu + v))in,r'-iAu)in,r''=j2i ] {Av){n, * iAu){n, -T^'-'^ 

1=1 



I 
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Using again Proposition 3.23| (viii) and Proposition \i.26 we obtain 



\Yiu,v)\\.,^ < r^E^E2'l|A.||Sl|A.||tJ||A 

k=2 1=1 



1 ^ IcK I 

< 1 — r V -r, 7TT (2max(||An||o,^, \\Av\\o,w)f^^ 2\\Av\\a,w 

l"i| ^ - 1)! 

< SCF.cmaxdlullctt,, ||f ||o,t«)||At;||^,^. 
Observing that ||Aw||(j,tii < (1 + e'^)||v||o-,M; the claim follows. 

(iii) The proof of differentiability for Fi (as well as F2 and F3) uses the fact that 
these functions are sums over / and k,l < k, of monomials of the form {Av){n, •)*' * 
(Aii)(n, •)*('^~'), Therefore it suffices to first prove the continuous differentiability of 
each term in the sum and to show secondly that the sum of the derivatives converges 
uniformly in the corresponding norm. 

Because of the simple algebraic rules for convolution (see Proposition 3.23| ) it is 
straightforward to check that for / > 1 the map 

Fi-.v^ {Av)in, * (A'u)(n, •)*^''"'^ 

is a map from Ca,w into Ca,w with derivative 

DF4{v)x = l{Av){n, O*^'"^^ * {Au){n, * (Ax)(n, •)• 

Proposition 3.23| (viii) yields the estimate in the corresponding operator norm 

\\DF^{v)\\ < (l + e'^)/(max(||Au||o,^,||At;||o,»))'=-' 



and with Proposition |3.26| we conclude the uniform convergence of the sum, as 



1 00 I \ ^ I h \ 

AE^E , (l + e'^)U2max(|H|o,.,||HM)^-^ 

fc=2 1=1 \ ^ J 

< ^ E 71^(1 + (4max(||n||o,., 

< 2(1 + e'^)C'i7,c4max(||n||o,«,, |b||o,-u;)- 

This proves everything about the first derivative. For the second derivative we can 
proceed similarily. We get 

D^F4{v)[x, y] = l{l - l){Av){n, * (An)(n, * (Ax)(n, •) * {Ay){n, ■). 
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The convergence of the sum is guaranteed by 



1 



k=2 1=1 



(1 + e")2/(/ - 1) (2max(||n||o,^, \\v\\o,^)) 



k-2 



-. oo 



(1 + e'')24(4max(||n||o,»,||u||o,^)) 



(fe-2)! 



k-2 



(iv) The proof is rather similar to the one just given. For A; > 2, / > 1 let 

-F5 is a map from £o,w) iiito Cfj^yj with 

DF^{u)x = (A;-/)(Au)*(*^-'~i) * (Au)*' * Ax. 
\\DF^{u)x\\„^yj < (k-l) (2max(||u||o,u,, ||f ||o,m«))''"^ (1 + e'^)lbll<7,M;2||x||o,«,. 



Proposition 3.26 gives the uniform convergence of the sum. The remaining part of 
(iv) can be easily seen from Proposition 3.2(:, Proposition 3.23| (viii) and the just 
given formula. 

(v) Applying the procedure again, we first convince ourselves that for I > 1 the 
function 

Fe-.v^ iAu)<^-^-^'^ * (Av)*' * Ax 
is a map from Ca,w into C(^-a-',w with derivative 

DFe{v)z = l{Au)<^-^-^^ * {Av)<^-^^ * Ax * Az. 

The sum of the operator norms of the derivatives can uniformly be estimated by 



-. 00 



2^^ (2max(||tx||o,^, Mo^n,)? ^ 2||x||_,/,^(l + e") 



\a,\f^^{k-2)\ 
< 8(l + e")C'F,c||x||_^/,^. 

This concludes the proof of the proposition. 



O 
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3.2.4 The equations for the Fourier coefRcients 



Following the ansatz described in Section |3.1| , we are now ready to give sufficient 
conditions for the Fourier coefficients u{n,m) and v{n,m) in order to obtain real 



solutions for the driven nonlinear lattice described by equations ( |3.1[ ) and (3.2). 
Recall from equation (^) the definition (5m = — 2 + -pvp^- 

Lemma 3.31 Let F,c,^, {bm)mGZ,'w satisfy the general assumptions. Suppose there 
exist u,v £ Cq^w, for which the following conditions hold. 

(1) 

Vn > 1, m G Z : u{n — 1, m) + 5mu{n., m) + u{n + 1, m) + W{u){n, m) = 0. 



Vn > 0, m G Z : u{n, —m) = u{n, m). 

II II / PP^^ 
o 

(2) 

Vn > 1, m G Z : v{n — 1, m) + 5mv{n, m) + v[n + 1, m) + W{u, v){n, m) = 0. 



Vn > 0, m G Z : v{n, —m) = v{n, m). 

II II / Z'^'': 
o 

(3) Vm / : e6m - t;(0, m) - u(0, m) = 0. 
Then the family of real valued and periodic functions 

Xn{t) := cn + ^ a(n, m)e^"^^\ for n > 1, 

a(n, m) := n(n, m) + v{n, m) + efem — ""(0, m) — v{0, m), for n > 0, 
solves the equations 

where 

xo{t) = e J2 bme'"^^'. 
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Note that the sequence a is also defined for n = and that a(0, m) = ebm, the 
Fourier coefficients of the driver. 

Proof : First we note that the twice differentiable functions for n > 1. 

In fact, we know that {5mu{n,m))mei, and {5mv{n,m))mez are sequences in £i , as 
they can be expressed in terms of the £i sequences n(n — 1, •), u(n + !,•), W{u){n, •) 
and v{n — l,-),v{n + 1, ■),W{u,v){n, ■) . But this imphes that {m?a{n,m))m^i^ is 
in ii, which yields the regularity of x„. Furthermore it is immediate that all 
functions x„ are real valued and periodic with period 

Let us now turn to the main point of the proof, namely to verify that the Xn 
are solutions of the driven lattice. As Xn and F{xn-i — Xn) — F{xn — Xn+i) are 
both continuous functions of the same period, it suffices to show that their Fourier 
coefficients coincide. One checks from the definitions that 



Vn > 1 : Xn-i{t) - Xn{t) = -C+J2 He*""^* 

= -c+ ^ A(u + t>)(n,m)e™^*. 



Substituting into the Taylor series for F yields 



F{Xn-l{t) - Xn{t)) = 

k>0 



fc>0 Vmez / 

EfrE(^(^ + ^))(^'-r'Me™^* 

fc>0 mez 

E (E^(^(^+^))(^'-r'Mje^'""*, 



mez \fc>0 



where all the manipulations are justified as the sums converge absolutely (compare 
with proof of Proposition 3.29| ). We can now read off the Fourier coefficients. 

^ [F{Xn-l{t)-Xn{t))-F{Xn{t)-Xn+ime~''^^'dt 

= ai [A(u + v)(n, m) — A{u + v){n + 1, m) + W{u){n, m) + W{u, v){n, m)] . 
On the other hand, using condition (3) of the hypothesis 

— [ " Xnit)e-'"'^'dt = -{-fmfa{n,m) 
2tt Jo 

= — (7m)^(n + f )(n, m). 
The equality of the Fourier coefficients follows from (1) and (2) of the hypothesis. 
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It is a simple corollary of the proof of the Lemma, to see that condition (1) is 
satisfied if we have a "small" solution of the doubly infinite lattice equation. This 
is stated more precisely in the following remark. 

Remark 3.32 Suppose that u : Z x Z ^ C satisfies Z^mez ^^Pnez \u{n,m)\ < ^-|^. 
For n £ Z set 

Then u satisfies condition 
(1') 

\/n £ Z,m £ Z : u{n — 1, m) + 5mu{n, m) + u{n + 1, m) + W{u){n, m) = 0, 



Vn G Z, m G Z : u{n, —m) = u{n, m), 
if and only if Xn^ (t) is a real valued solution of the differential equation 

Xn{t) = F{Xn-l{t) - Xn{t)) - F{Xn{t) - X„+i(t)), for U £ Z. 



3.3 Solving for the non resonant modes 

The present section is devoted to the proof of our basic result, which was explained 
and motivated in Chapter |l| and in the introduction to this chapter. Before the 
theorem can be stated we recall some notation of Section |3.1| and we add a few 
definitions. 

In equation (^^) we have set 

^- = -^ + F^y 

Furthermore we denoted in equation (|3.8|), 



mo = max{m £ Z : < < 4}. 

By separation of variables one obtaines solutions of the free linearized problem of 
the form 

= m£Z, 
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where 

(3.33) zl^ + 5mZm + 1 = 0. 

The case \m\ > ttiq corresponds to 5m > 2 and therefore we can pick z,m to be the 
sohition of the above equation with \zm\ < 1, which is given by 



(3.34) V|m| > mo : := + y ^ - 1, for J„ > 2. 

In the case < |m| < mo, the general assumptions on the frequency 7 imply 
that \5m\ < 2. We choose for Zm the solution of equation ( p. 33 ), which corresponds 
to an outgoing wave yn(t), i.e. Zm = e*^™ with 



(3.35) Pm = — sgn (m) arccos 



' 2 

This explains the choice of the sign in equation ( |3.9[ ). Note, that for all < |m| < tuq 
we have 

(3.36) |sin/3^|>0, 
as 1 5m I < 2. 

Definition 3.37 Let 7,F, c, {bm)m£z,w satisfy the general assumptions. 

'^ln(|^mo+l|))- 





:= min( 




1 




:= ^ao. 




:= max 



+ 



1 - e-^^'o \(>(To _i 1 _ g-2<To 
1 

max < — — - — - : < m < mo 



1 - e-'^i 1 1 sin P, 
Pu ■= w{m)\u{0,m)\. 

\m\>mo 

Remarks : 

• The constant Ck is well defined, as uo, cri > and \sin(3m\ > for < |m| < 
mo. We will see below, that Ck is an upper bound on a linear operator K 
which is related to the to the inverse of the operator 0mez -^m (see Proposition 



3.50| below). 
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• The use of pu will not become clear before Chapter ^. There we will see that 
this quantity is the main tool for proving that the sequence v ( in the notation 
of Lemma 3.31| ), which will be constructed in the next theorem, is of higher 
order. 

The following is the main result of this chapter. 

Theorem 3.38 Let 'j,F,c,{bm)m.GZ,'w satisfy the general assumptions. Further- 
more we assume that there exists a choice of constants G N, cq > 0, Cq > 1, and 
a map {q G : \q\ < cq} jCq.w ■ Q ^ u{q), such that 

(1) M\q\ < Co : ||n(5)||o,^ < ci and Pu{q) < ^, where 

1 PF,c\ 



(3.39) 



ci := mm 



(2) q 1-^ u{q) is a map from {q G : \q\ < cq} to £-ai,w and the following 
estimates hold. 



y\q\<co,l<j <N : 



_d_ 

dqj 



u{q) 



<Co. 



-cri,w 



(3) For all \q\ < cq, n > and m G Z : u{q){n, —m) = u{q){n, m). 

Then for all {q,e) G M^+^ with \q\, \e\ < min(co, ■^^) there exists a unique v G /2o-o,«) 
with the following properties (i)-(iii). 

(i) Mm G Z, n > 1 : v{n — 1, m) + 6mv{n, m) + v{n + 1, m) + W{u{q),v){n, m) = 0. 

(ii) \\v\\ao,w < HPuiq) + |e|) < ci < 
(Hi) V|m| > niQ : v{0, m) = ebm — u{q){0, m). 
Furthermore the following holds. 



(iv) For all n > and m G Z : v{n, —m) = v{n, m). 
(v) The map {q,e) i — > v is a map into La^^w 
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Proof : The proof proceeds via a Banach fixed-point argument for v and 
the derivatives of v with respect to the parameters qj,€. First we define a map T 
(Step 1), which we then show to be a contraction on a certain set (Step 2). We 
conclude that the first component of the fixed-point of this map is a solution of the 
equations given in (i) and (iii) (Step 3). Step 4 settles the question of differentiabil- 
ity for V and Step 5 deals with the remaining properties (ii) and (iv). 



Step 1: Definition of the contraction. 



First we will turn equation (2) of Lemma 3.31 into a fixed-point equation v = 
T{q, e, v), depending on the parameters q and e, by applying the inverse of 0mez 
on it. As it was pointed out in the introduction of the present chapter, the oper- 
ators Lm are not invertible for \m\ < tuq. Nevertheless we will define a formal 
inverse for \m\ < niQ, acting on exponentially decaying sequences, as motivated 
in the introduction. Note the special role of m = 0, where the Green's func- 
tion of Lq grows linearly. Using the fact that the nonlinear term is given by 
W{u,v){n,0) = Y{u,v){n,0) — Y{u,v){n + 1,0), we end up with a bounded ker- 
nel acting on Y{u,v){n,0). We then proceed to define the maps Tq^^j, which give 
rise to the fixed-point equation for v in their first argument and to the fixed-point 
equation for the partial derivative of v with respect to the j-th component of the 
parameters in their second argument. The map T is introduced in order to show 
that the fixed-point v of the map T depends smoothly on the parameters q,e in 
Step 4. 

Let u := (g, e) denote the parameters in the construction. Furthermore it is 
convenient to scale these parameters by a factor 

(3.40) 7/:=min(co, 



4Co 



Hence we can choose v £ U, with 

(3.41) U := {{q,e) G R^+^ : \q\, \e\ < 1}. 
Let us further define two sets on which the map will act. 

(3.42) B := {v £ C^^^w ■ \\v\\ao,w < ci}. 

(3.43) B' := {y G Cao-2ai,w ■ ||y||(7o-2cri,«; < Cl}. 

We now define the map T^v, v) explicitly. For v £ U,v £ B \et 



48 



• m = 0,n > : 



T{v,v){n,Q) := - ^ Y {u{riq),v)(k,^). 
k=n+l 



• < \m\ < mo, n > : 



T{i',v){n,m) := — — — ^ W{u{r)q),v){k,m) sm{n — k)P„ 



• \m\ > mo, n = : 

r(i^,v)(0,m) := 77e6^ - u{riq){0,m). 

• \m\ > mo, n > 1 : 

oo 1 2min(n,A;) 

r(i.,i;)(n,m) := £ ^T^^ zj^J-^^l+^W^Kw), m) 

J- 



k=l 



+ 2;;;,(?7efe^ - t/(77g)(0,m)). 

It is useful to rewrite T in the following way. Define the linear map K, which acts 
on spaces £1^,^. 

oo 

(3.44) (Ky){n,m) := K{k, n, m)y{k, m), 

k=l 

with kernel K{k, n, m). 
• m = : 

K{k,n,0) := 



0, for k <n 
— 1, for k >n + l. 



• < |m| < mo : 

K{k, n, m) := 

• \m\ > mo : 



0, for k <n 

^ [sin(n - k)Prn — sm{n + 1 - A;)/?^] , for A; > n. 



I 2min(n,fe) -. 2min(n,fc— 1) 

K{k,n,m) := ^ ' '"^^ , 4-^1+^ - ' ~ , 
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It is straightforward to check that 

(3.45) T{u,v){n,m) = [KY{u{r]q),v)]{n,m) + {ijebm - u{riq){0,m))zl^l{\^\y^^y. 

Now we can define the map T, which acts on the complete metric space B x B', 
equipped with the norm 

(3.46) \\{v,y)\\BxB' := max(||t;||^(,^^, ||y||^(,_2CTi) • 

Fixu eU,l < j < N + 1. 
fuj : Bx B' — > Bx B', 



(3.47) Z 



^,3 



V \ I T{v,v) 
y j ' y hj{u,v) + K{D^Y){uir]q),v)y 

where we now give explicit expressions for hj{i',v). There are two cases. 
Case 1: i^j = e 

(3.48) hj{u,v){n,m) := r]bjnzl!a'^{\ra\>mo}- 
Case 2: i^j = qj. 

(3.49) hj{u,v){n,m) := ri[K{DuY){Dju)]{n,m) - r]{Dju){0,m)zl^l{\^\^„^^y 
This completes the definition of T. 

Step 2: Tjyj : B x B' — > B x B' is a contraction. 

We first obtain a bound on the norm of the linear operator K. 

Proposition 3.50 For all ai < a < ao, the linear operator K maps Ca^w into La^w 
and the corresponding operator norms of K are hounded by Ck- (See Definition 



Proof : The proof is a consequence of the following estimates. 



• m = 0, n > 



|(i^y)(n,0)|e'^" < e'^" e"'^' sup |e'^^y(j, ( 

k=n+l 

< , '^■^ sup|e"^y(j,0)|. 
1 - e '^i j>o 
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< \m\ < rriQ, n > : 

\iKy)in,m)\e" 



< 



< 



J2 e--'=sup|e--'y(j 

k=n+l 

sup |e y{j, m 



1 - e-"^ sin /3, 



m j>0 



\m\ > 1710,12 = : K{k,0,m) = 0. 

|m| > mo, re > 1 : By definition of uq, \zm\ < e~^°"° (see Definition 3.37 ), 
whicli yields the estimate for the Greens function 



< 



_£__^_g 2(To|n Convolution of this bound with an exponentially decay- 
ing sequence gives 

oo n— 1 oo 



-2ao(n—k)—(jk 



+ E 



fc=i 



fc=0 fc=n 
g{2(To-cr)n _ 2 



2(Ton 



g(2o-o-o-) _ I 



+ e 



-2o-o{A;— n)— crfe 



1 



1 



+ 



go-o _ 1 1 _ e-2ao 



and we arrive at 

\{Ky){n,m)\e' 



1 _|_ e~2'^o /I 1 \ 
< , ^ T + -. ^ sup|e"^y(j,m)|. 



O 



Now we are ready to prove that Tj^j- is a contraction. We use the estimates of 
Proposition 3.29| and equations ( 3.45| )-( |3^ ). 



2 , Ci Ci 3 



\\T{u,v)\U,^ < CkCfA + J + J < 4C1. 



Furthermore it was shown in Proposition 3.29| (iii) that, for the u, v in question, 
the map v 1— > Y{u, v) is C"^ from Caa^w into Cuq^w and the ususal operator norm 
of the derivative is bounded by Cp^cCi- Therefore 



\\T{u,v')-T{u,v) 



\ao,w < CkCf,cCi\\v' - v\\ao^w < ^1^' " v\\mhW- 
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We have for u e U,v e B,y e B', that 

(JQ —1(7\ ,W 

The bound on the second derivative of Y with respect to v (see Proposition 



3 .291 (iii)) allows us to make the following estimate. 



\\K{D,Y){u,v')y' - K{D,Y){u,v)y\U.2a^,v, 

< \\K[iD,Y)iu,v') - {D.,Y)iu,v)]y'\U-2a^,w 
+ \\K{D,Y){u,v){y' -y)\\ 

< CkCf,cCi\\v' - ■u||cro,«> + CKCF,cCl\\y' - y||ao-2fTi,ii) 

< ^ max {\\v' - w||<xo,t«, \\y' - y|Uo-2ai,t«) • 

Finally we are dealing with the estimates for hj{v,v). Corresponding to the 
definition (|3l^ ), ( p9| ) we have to distinguish two cases. 



Case 1: i^j = e 



\\vbmZm'^{\m\>mo}\\o-Q-2cn,w < ^ < "tCi. 



4 

Case 2: I'j = qj- 

Using hypothesis (2) of Theorem 3.38| we see immediately that 



II^Pi^^)(0,H2;ml{|m|>mo}IUo-2<Ti,«; < V^O < \cx. 



From Proposition 3.29| (iv), we obtain 



\\r]K{D^Y){u,v){Pju)\\„,^2a,,^ < 2r?C/rC7^,eCiCo < ^ci. 

o 

Finally we use that for fixed x G C-ai,w, the map v ^ {DuY){u, v)x is from 
Ca-Q^w into Co-Q-2c7i,w with a bound on the derivative as given in Proposition 
g (v). 



\\r]K[iDuY)iu,v')-{D^Y)iu,v)]iDju)\\„,^2aun. < t^CkCf^cCoWv' - v\ 



The claim of the present step is a consequence of all these estimates. Thus we 
have proven the existence of a fixed-point of T,yj in B x B'. 
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Step 3 : For v £ B the following eqivalence holds. 
V satisfies properties (i) and (iii) of Theorem 3.35 <^=^ T{h', v) = v. 



Recall that we have scaled by the factor r] in the beginning of the proof. Then 
property (iii) is evident as T(z^, v) satisfies it by definition. Therefore it suffices to 
prove that for al\ v £ B,n > l^m £ Z: 

T{i', v){n — 1, m) + 5mT{v, v){n, in) + T(i/, v){n + 1, m) = —W{u, v){n, m). 

We will verify this by evaluating the lefthandside for all different cases. 

• m = 0, n > 1 : 

oo CO oo 

LHS = -Y,Yiu,v)ik,0)+2 J2 Y{u,v)ik,0)- ^ Yiu,v)ik,0) 

k=n fc=n+l k=n+2 

= -Y{u,v){n,Q) +Y{u,v){n + l,id). 



< |m| < mo, n > 1 : 



LHS = — -7^ ( W{u,v)ik,m)Gik,n,m) +Wiu,v){n,m)s\iii—(3ri 



vn 



where by the definition of /3m (see ( ^.9D ) 

G{k, n, m) := sin(n — 1 — k)(3m — 2 cos Pm sin(n — k)Prn + sin(n + 1 — k)Prn- 

As sin(n — 1 — k)Prn + sin(n + 1 — k)Pm = 2 cos Pm sin(n — k)fjm, we conclude 
G(/c, n, m) = 0. 

\m\ > rriQ, n = 1 : 

LHS = {l + 6mZm + z'^){riebm-u{r]q){0,m)) 

00 / 2min(2,fc) \ 

+ J2 W{u,v)ik,m) 5mzi + ~ , • 

k=l \ ^ J 
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The identity 1 + SmZm + Zm = ^ implies that all terms in the sum vanish with 
the exception of A; = 1. In fact 

• k = l: 

, 2min(2,A;) 
S I 2|2-fe|+l _ X . 



• k>2: 



, 2min(2,fc) 
J- 



• |m| > mo, n > 2 : 

LHS = (zl-' + 5mzl + zZ+^)iriebm-u{riq){0,m)) 

oo 

+ ^ m)G{k, n, m), 

k=l 



with 

G{k, n, m) 



-. 2 min(n— l,fc) 2min(n,fe) 

^~ Zm e ^ — Jn-k\+l 

1 _ ~2 ^rn 1" 1 _ ^2 

I 2 min(n+l,A;) 

I -*- ~ ^"-^ |n+l-fe|+l 

1 _ ~2 

J- -O™ 



We show that G(A;, n, m) = —5k,n- To that end we evaluate (1— z^)G(A;, n, m) = 

• for < n: 

= (1 — Z^)zl^ '^(1 + 5m^m + ^m) ~ ^■ 

• for k = n: 

• for k > n: 



Let V E B satisfy (i) and (iii). Denote d := v — T{v, v). We want to show that d = 0. 
The following list of properties for d is immediate. 
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(I) d e Cao,w 

(II) d{0,m) = for \m\ > rriQ (as v satisfies (iii)). 
(Ill) Vm G Z, n > 1 : d(ra - 1, m) + 5md{n, m) + d{n + 1, m) = 0. 

The last relation enables us to express d{n, m) in terms of (i(0, m) and d{l,m). 

• m = : d{n, 0) = nd{l, 0) — (n — l)(i(0, 0) . As d has to decay exponentially 
with n — ^ GO there is no other choice than Vra > : d{n, 0) = 0. 

• < \m\ < mo : 

sinn/?^ sin(n-l)/3„i 

d[n,m) = — : — - — a(l,m) : — a(0, m) 

siufim sin/?^ 

_ — n^m (^^(^^^ jy^-^ _ ^(^Q^ yji^ (,Qg _j_ cosn/3TOC^(0, m). 
sin/3^ 

/?m € (0, tt) and again the exponential decay of d(-, m) force d{Q, m) and d(l, m) 
to be zero and hence d{n, m) = for all n > 0. 

• \m\ > rriQ : Recall that we know already by (II) that d(0,m) = 0. Then 

gain we conclude that d(l,m) = 0. 

This concludes the proof of Step 3. 

Step 4: Differentiability of v with respect to the parameters u. 

It is a standard and well known problem to prove smooth dependence of the solution 
of a contraction problem on the parameters. There is a small technical problem as 
ti is a C^-function of the parameter q only in the space and not in the space 

JCq^w But, at least formally, differentiating v = T(i',v) with respect to Vj gives 

^ = {l-D^T)-^^T{u,v). 

The lack of differentiability of u{q) in Co,w prevents that -^T lies in CaQ^w , but 
it lies in CaQ-2ai,w Fortunately, however, DyT maps any Ca,wj cri < a < ctq, into 
itself with small norm. Hence v is differentiable in an appropriate norm. We find it 
convenient to proceed as follows. 
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The fixed-point {v,y){v) of T^j can be constructed as the hmit of the iteratives 
of the map, i.e. let 

{v^{v),y^{v)) := (0,0) 

and define inductively 

{vs+i{v),ys+i{v)) ■■= f^j{vs{v),ys{v)), 

then {v{u)^y[v)) = lims_>oo(us(z^), ys(z^)). The limit is in the norm of B x B' and 
uniform in v. 

We will prove inductively for each choice of z^j and for all s > that the following 
properties hold. 

(I) U — > £cto-o-i,w) : ^ Vsii') is continuous. 

(II) For each variable i^j, the map i^j Vg^i^) is differentiable as a map into 
■^CTo— 2(71,10 and 

-rT-Vsiy) = ys{v)- 

OVj 

(III) U — > CcrQ-Zcn,w '■ ^ ^ Usii^) is continuous. 

Once we have established (I)- (III), the proof of the claim is immediate. In fact, we 
can deduce that for all s G No : ^ 'Us(z^) is a map from U into Cao-3ai,w = ^cti,w 



(see Definition 3.37 ), with partial derivatives = ys- The uniform convergence 



of {vs,ys){'^) to the fixed-point of T^j, (t;,y)(i/), in the norm oi B x B' yields the 
desired information, that u i— > v{i') is a map from U into Cai,w and that ^ = y. 

Let us therefore return to the statements (I)-(III). They are trivially satisfied for 
■s = 0. We will now prove the induction step s — > s -|- 1. 

(I) Continuity: 

Let u, v' S U . We have to show that 

(3.51) \T{y' ,Vs{v)) -T{y,Vs{y))\aQ-ar,w — > 0, cSSV ^v. 

To simplify the notation, denote v := Vs{i'),v' := Vs{i''),u := u{riq),u' := u{riq'). It 
is not hard to verify that 

\\[{r]ebm-u {0,m)) - {rjebrn-u{0,m))]zl^l{\rn\>mo}\\'7o-(7i,w < vie -€\ + \\u -u\\-ai^w 
(3.52) 
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Expressing the difference by telescoping sums, one obtains for I > 1 

||(An')*('-') * {AvY - (A^.)*(*^-') * [AvY'lU^^,,^ 
< (1 + e"«)A:(2ci)^-^ {\\u' - u\\-^,,^ + \\v' - v\\ao-a„w) • 
Proposition |3.26| and Proposition |3.5[l| imply 

\\KY{u',v') - KY{u,v)\U.^,,^ 
(3.53) < 8(1 + e''°)CKCF,cCi {\\u' -u\\^^,,^ + \\v' -v\\^,-^,,^) . 

Equations ( ^.521) , (|3.53| ), assumption (2) of the Theorem and the induction hypoth- 
esis suffice to prove equation ( p. 51 ). 



(II) Existence of partial derivatives: 

We consider only the more difficult case vj = qj. The proof for the case I'j = e 
requires only a proper subset of the arguments given below. Denote ly' := {q + hqj, e) 
and let u', w', u, f , Dju' , DjU, Djv' , DjV have the obvious meaning. We have to prove 



(3.54) 



lim 



-{T{u',v')-Tiu,v))-ys+i{u) 



= 0. 



We break this statement up into several estimates. 

— (n'(0, m) — u{0, m)) — r]{Dju){0, m) 



fjQ — 2cri,U' 



^m-'-{|m[>mo} 



cro — 2cri 



(3.55) < 



h 



[u' — u) — rjDjU 



which tends to as /i ^ by assumption (2). Next we look at the monomials, from 
which Y, DuY and DyY are built up. We use the induction hypothesis, which says 
that Us = DjVg. 

^ [{Au')*'^''-^^ * (At-')*' - (An)*^^-') * {Av)*^' 
-{k - /)(An)*('=-'-i) * (Av)*^ * r]A{Dju) 
-Z(An)*('=-') * (Ar;)*('-i) * A{Djv) 
= i ((An')*^^-'^ - (Au)*(^-')) -{k- 0(An)*('=-'-i) * r]A{Dju) * (At- 

+{k - /)(An)*('=-'-i) * T]A{Dju) * (^{AvY - (Av)*') 



r\*l 



+(An)*('^-') * 
(a) + (6) + (c) 



((At;')*' - (At;)*') - l{Av)<^-^^ * A{Djv] 
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The three terms are now investigated seperately. Assume I > 1. 
(a) Telescoping differences twice we obtain 

i ((An')*^^"'^ - (Au)*(^-')) -{k- /)(An)*('^-'-i) * r]A{Dju) 



h 

k-l-lj-l 

J2 II(An')*^ * (An)*('=-'-^-2) * A(n' - u) * A 

j=l i=0 



u — u 



+ {k - l){Auf^-^-^^ * A ( ^^-^ - rjDjU 



This imphes, that 



(a)|Uo-2<xi,«; < 2(l + e'^«)A;(A:-l)(2ci) 



fc-2 



X m — u\ 



u — u 



h 



+ Cl 



— CTl ,10 



u — u 



h 



rjDjU 



-a\ ,w , 



(b) Recalling the definition of Cq in the statement of the theorem, it is easy to see 
that 



CTo-2<xi,«, < 2(1 + e'^°)A:(A: - l){2cif-^iiCo\\v' - v 



(c) Proceeding as in (a) we obtain 

II \pJ l|ci'o~2(Ti,to 

< 2(l + e'"o)A:(/c-l)(2ci) 



fc-2 



X Ih - v\\o^^ 



V — V 



+ Cl 



o-0-2(Ti,u) 



V — V 



DjV 



o-0-2(Ti,«) , 



Substituting all these estimates in the power series and using Proposition 3. 26, the 



induction hypothesis, assumption (2) and equation (3.55) we arrive at the assertion 
of equation ( |3.54D . 



(Ill) Continuity of partial derivatives: 

Using the notation and the methods of the last two proofs, we see that the fol- 
lowing estimates are enough to prove the claim. 



||r?((I),V)(0,m)-(D,n)(0,m))z-l 

{|m|>mo} lUo— 3' 

< r]\\DjU - Dju\\^^^^u}- 
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CTQ — 3cri,ui 



The usual telescoping technique readily yields for I > 1 the following estimate. 

{Au')^''-^-'^^ * (At;')*' * A{Dju') - (Au)*^''-^-'^^ * (Av)*' * A{Dju) 
< 2(l + e'"«)^(2ci)''^^ 

X {Co\\u' - + Co\\v' - v\\ao~auw + Cl\\Dju' - DjU\\^ai,w) ■ 

• We know from Step 2 that ||-C'jt)'||o.f,_2(Ti,«), ll-^^i^llo-o-2(ri,«; < ci. Proceeding as 
above we obtain 



cro~3ai,w 



{Au')*^''-^^ * (Av')*^'"^^ * A{Djv') - (Au)*(^-') * {Av)<^-^^ * ^{Djv) 
< (l + e"")A;(2ci)'=-^ 

X {\\U - u\\-ai,w + \\v' - v\\ao-auW + \\Djv' - DjV\\ao-3ai,w) ■ 

The proof of Step 4 is completed. 

Step 5: The remaining properties (ii) and (iv). 



(ii) is a consequence of the fact that for all G C/ the map v i— > T(i/, v) sends 
{v G C^o^w '■ \\v\\ao,w < 2(r/|e| + Pu{riq))} into itself. In fact, from Proposition [3.25 
and equation ( 3.45| ) it follows, that 



\\T{u,v)\\ao,w < CKCF,cCl\\v\\ao, w + ivl^^l + Pu{rjq)) 
< + Pu{m))i'^CKCF,cCl + 1). 

Note that the appearance of rj in the proof, which is not present in the formulation 
of property (ii) in the theorem, comes from the scaling of the parameters which we 
performed at the beginning of the proof. 

(iv) can be shown inductively, following the iterative construction of the fixed- 
point which we have already employed in Step 4. vq = clearly satisfies the reality 



condition, i.e. Vo{n,—m) = vo{n,m), for all n > 0, m € Z. Using Propositon 3.23 
(v) it follows that Vs+i satisfies the reality condition as {bm)mez,u{r]q) and Vs do. 
This property is preserved as we pass to the limit s — > oo. 

O 
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3.4 An immediate application: high frequency driving 



Suppose that mo = 0, or equivalently 7^ > 4F'(— c). In this case we can choose u{q) 
identically equal to zero in Theorem ^.38 and we obtain v{e) with the corresponding 
properties. If we now take these choices for u and v we see immediately that all the 



conditions in Lemma 3.31 are satisfied and hence we have constructed a periodic 
solution of the differential equation given by ( |3.1| ) , (^) . This proves the following 
result. 

Theorem 3.56 Let F,c,{bm)m&z,'w satisfy the general assumptions. Furthermore 
let 7^ > 4F'(— c). Then there exists a neighborhood D of 0, such that for all e ^ D 



there exists a sequence v{e) G £0-0 ,«> with ao > as defined in Definition \3.37 such 
that 

Xn{t) := cn + e6o - w(e)(0,0) + ^ i;(e)(n, m)e^"''^*, for n > 1, 



is a time periodic, real valued solution of the the differential equation, given by (3.1) 



and ( \3.^ ) 



Remark 3.57 We have restricted our attention to the case where F'{—c) > 0. // 
F'{—c) is negative, however, then 7^ > AF'{—c) for all 7 G and it is easy to see 



that Theorem 3.56 holds without any restrictions on the driving frequency. 
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Chapter 4 



General Lattices 



In the last section of the last chapter we have seen that in the case of uiq = 0, it is 
possible to construct periodic solutions for arbitrary lattices, where the interacting 
forces between neighboring particles satisfy the general assumptions. The goal of 
this chapter is to obtain the same result in the case of mo = 1, i.e. 

(4.1) F'(-c) < 7^ < 4F'(-c). 

The difference from the case mo = is that now Theorem 3.38| no longer yields a 
sequence v which solves ebm = m(0, m) + f (0, m), for all m ^ 0, but only for |m| > 1. 
Therefore we have a resonance equation for m = 1. In order to be able to solve this 
complex valued equation ( for m = — 1 the equation is the complex conjugate of 
the equation for m = 1), we must obtain a family of sequences u{q) depending on 
two real parameters qi and q2- The idea is to obtain u(q) by constructing travelling 
wave solutions for the doubly infinite lattice. Physically one may view these as the 
waves the driver excites and which travel through the lattice. They can be observed 
almost unperturbed away from the boundary some time after we let the driver act 
on the system (see Figures C/7 and C.lOl ). 
More precisely we make the ansatz 

(4.2) x(f ) (t) := cn + ^ ^^^ym{i3n+^t) _ 



In Section |4.l| we will give an equation for r which we then solve by a Lyapunov- 
Schmidt decomposition. The idea is as follows. It is easy to see that the linearized 
equation at r = is given by a diagonal operator A(/5) with entries 

(4.3) A(/3)(m,m) := A(/9,m) := 2 cos /3m + 6m- 
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If we choose j3 = (5i (see ( ^.91 ), /3i exists as \5i\ < 2) , then A(/3) has a nontrivial 
kernel and we can apply the decomposition procedure. The assumptions on 7 as 
given in ( [4.1D imply that 5m > 2 for \m\ > 1 and therefore A(/3)(m,m) is bounded 
away from zero for arbitrary /? G M. Hence the infinite dimensional part of the 
decomposition poses no problems. The finite dimensional part needs additional 
consideration. First, one has to use various symmetries to obtain the correct count 
of variables. Then by expanding the degenerate equation to second order one proves 
that the finite dimensional part can be solved by choosing the spatial frequency /3 
as a function of the two remaining real parameters qi and (72 with 

\Piq)-(3i\<C\q\^ 



In Section O we then take these solutions r{q) and show that 



(4.4) n(g)(n,m) :=r(g)(m)e^'^'^^('?) 



satisfies all the conditions of Theorem 3.38| in Chapter ^ Note, that from equation 



( |4.4| ) it is already obvious that the map q^ u{q) cannot even be continuous in Cq^w, 
if (3{q) is not identically constant. This is why we had to introduce the weaker space 
C-a-^^w in the previous chapter, which led to some additional difficulties. Further- 
more, as we will see below, it will be necessary to relax the weight function w hy a. 
factor Tin?. We will anticipate this loss and construct r(m) in ii^w-, with 

(4.5) w{m) \= w{m){l + \m\f' . 
Note that w is again an admissible weight function. 

4.1 Construction of travelling waves 

We make the following definitions. 

(4.6) [A(/3)r](m) := A{l3,m)r{m), with A(/?,m) = 2cos(/3m) + 5m- 

(4.7) Ar(m) := (e"^'^"' - l)r(m). 

CXD 

(4.8) W(r)(m) := (1 - e'^"")— V ^(Ar)*^(m). 
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Finally, define an operator acting on ii^u) in the following way. 
(4.9) mr){m) ■.= e'^"'r{m). 

This multiplication operator on a Fourier sequence corresponds to a translation 
of the function, an operation under which the autonomous differential equation is 
invariant. 



The next lemma follows easily from Remark 3.32. 



Lemma 4.10 Suppose that /? E M and r : Z ^ C satisfies \\r\\£-^ < For n £ Z 

set 

Then r satisfies conditions (A) and(B), with 
(A) K{^)r + W{r) = 0, 



(B) Mm € Z : r{—m) = r{m), 
if and only if Xn'^ (t) is a real valued solution of the differential equation 

X'nit) = F{Xn-l{t) - Xnit)) - F{Xn{t) - X„+l(t)), for n £ Z. 

Before starting the construction of solutions of equation (A), we have to inves- 
tigate the smoothness and symmetry properties of the nonlinearity W. 

Proposition 4.11 There exists a constant Cf,c (which can he taken to he the same 
as in Proposition 3.2^) such that for all r G with \\r\\£^^ < the following 



holds. The series in the definition of W{r) converges absolutely to an element in 
ii^w- Furthermore 

(i) W : ^r £ £i^yj : < ^-|^| — > (i,w is a smooth map and the following 

estimates hold. 

\\W{r)\U^^,<CFjr\\l^. 
yxeh,n, ■■ \\DW{r)x\U^^^<CF,c\\r\\e,J\x\U^^^. 
Vxi,X2G^i,tD : \\D'^W{r)[xi,X2\\\ei^ < CF,c\\xi\\eiJ\x2\\ei^. 



(a) If for all m E Z, r(— m) = r{m), then we have W{r){—m) = W{r){m), for all 
m G Z. 
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(Hi) Wmr) = T^W{r). 

(iv) r{m) G iM for all m G Z implies W{r){m) £ iR for all m £ Z. 

Remark 4.12 Property (Hi) reflects that the underlying equation is autonomous. 
Property (iv) corresponds in the original space to the fact that if Xn{t) is a solution, 
then —X-n{—t) is again a solution. 



Proof : The methods for proving (i) form a proper subset of what has already 



been done in Section 3.2.3, Proposition |3.29, and therefore these arguments are not 



repeated here. The properties (ii) and (iii) can be deduced from the corresponding 
properties of the convolution (see Proposition |3.23| ). Finally we observe that r(rn) G 
zM implies (Ar)(m) E e^^^R. Arguing as in (iii) we conclude for all k > 1, that 



(Ar)*'^(m) G e~~™'R, which yields (iv). 



O 



For the construction of the solution we will use a Lyapunov- Schmidt decompo- 
sition. We introduce the projections 



r{m), for \m\ > 1 
0, else 



(4.13) r ^ (Pr)(m) := 



(4.14) Q:=I-P. 
Denoting cp := Qr, ^ := Pr, we obtain two equations: 

(4.15) A{P)fi + PWiip + fi) = 

(4.16) A{p)ip + QW{ip + fi) = 

For given ip and /? we will be able to solve the first, infinite dimensional equation 
for fx (see Lemma 4.19| below). Substituting this solution into the second equation 



and using various symmetries we end up with one equation which we can solve by 
choosing /? as a function of ip, with /? close to Pi and (p close to 0. As mentioned 
in the introduction of the present chapter the spatial frequency /?i was chosen, be- 
cause A(/?i, 1) = 0. It is obviously true that A(/3_i, 1) = (recall from ( |3.9[ ), that 
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Pi = — and therefore one could replace (3i by in the construction which 
follows. However, in all the numerical experiments that we have performed (see e.g. 
Figures C.7 and C.IO) we only observe solutions corresponding to This is related 



to the direction in which the energy flows in the lattice and is explained in more 



detail in Remark 4.23 at the end of this section. 

Let us now parameterize ip. Observe that the second equation ( [4.16| ) is automat- 
ically satisfied for m = (use (p!6|) and (|48|) ). On the other hand (Ar)(0) = 0, 
independent of the value of r(0), i.e. the value of r(0) has no influence and we can 
normalize it to be zero. This freedom reflects an additional symmetry, namely if 
Xn\t) is a solution of the differential equation of the doubly infinite lattice, then 
Xn\t) + constant is again a solution for any constant. Consequently we parame- 
terize only a subspace of the range of Q, namely 

^{qi + iq2), for m = 1 

(4.17) ip{qi,q2){m) := ^ }.(^q^-iq^)^ for m = -1 

0, else 

Under a slight abuse of notation we define : ^ R^, 



(4.18) 

We then have 




cos ^ — sin 
sin ^ cos 




Let us now solve the first equation ( |4.15| ) . 

Lemma 4.19 There exists a neighborhood U C R'^ o/(/3i,0,0), a map ji : U 
^i,uj, {P,Qi,Q2) "-^ /"(/S)?) CLiT'd a constant C, such that 

(i) A{^3)^liP, q) + PW{^{q) + /.(/?, q)) = 0. 

(ii) Qi2{p,q)=0. 

(Hi) V/3 : n{P,0) = 0. 

H MP,q)\U^^^<C\q\^. 

(v) y{q)+^ip,q)\\,^^^<!^. 
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(vi) iJ,{f3,q){—m) = fi{(3,q){m), for all m S Z. 
(vii) n{P,T^q) = T(.^x{(),q) . 
(via) ^{(3,0, q2) G «M. 

Proof : We show first that the existence of the function /i can be deduced from 
the implicit function theorem apphed to the map 

G:UxS^fie Ran (P) : ||^||,,^ < ^| ^ Ran (P), 

{P,q,^x) ^ fx + A{P)-'PW{v{q)+fx). 

Indeed, note that A{(3)~^ is well defined on Ran(P), as A(/5, m) 7^ for \m\ > 1 (fol- 



lows from condition (4.1) in the introduction of the present chapter). Furthermore 
it is easy to check that is a function of /5 in the operator norm, which shows 
together with Proposition 4.11 that G is a map, if U is chosen as a suitably small 



neighborhood of (/3i,0) in R . Using Proposition 4.11 again we see, that 

G(/3i,0,0) =0, 

and 

D,GiPuO,0)= idRan(P)- 

which allows us to apply the implicit function theorem and to obtain properties (i) 
and (ii) immediately, (iii) follows from the uniqueness of the implicit function the- 
orem and the observation, that V/? : G(/3, 0, 0) = 0. In order to show the remaining 
properties we will investigate the family of maps 

Note that a fixed-point of -f^/3,g corresponds to a zero of G. That -f^/3,g is a contraction 
for suitably small values of q and ^ is a consequence of Proposition 4.11| . In order to 



prove (iv) it suffices to show that Hp^g maps |/U G Ran (P) : WfJ-We^ ^ < C|(7p| into 
itself for a constant C, which will be chosen below. In order to verify this claim, let 
us denote 

1 1 

Ca := sup 



|m|>l 



-2 + Sm —2 + 62 
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Proposition 4.11 implies 



< C\q\^ 

if we choose 

C := ACaCf,c 
and we only allow g G such that 



'w{l)+w{-l] 



|2 -L ||/,||2 



vj{l)+w{-l] 



2CACF,cC\q\^ < 2- 

Property (v) follows trivially from (iv) by making \q\ sufficiently small. 
The last three properties can be proved in the following way. q) is the fixed- 
point of Hp g and can be constructed as the limit of the sequence ^kiP, q) as k ^ oo. 
We define inductively 



0, 

Hf3,q{fik{P,q)). 



(vi),(vii),(viii) are trivially satisfied for fiQ as well as for the sequence (p{q). Proposi- 
tion |4.11| permits us to conclude inductively that all fij^ possess the three properties, 
which then is preserved under the limit k ^ oo. 



O 



We turn now to the finite dimensional equation ( [4.16| ). 

A{PMq)+QW{^{q)+fiiP,q))=0. 

Let us first reduce the dimensions of these equations by factoring out all the sym- 
metries. We recall that the equation is satisfied for m = by (4.6) and ([4.^). Fur- 



thermore it suffices to solve equation ( 4.16 ) for m = 1 as the equation for m = — 1 
only yields the complex conjugate of the same equation. This leaves us with one 
complex valued or equivalently with two real valued equations. Finally we can make 
one further reduction using the invariance. It implies that a choice of {(3,qi,q2) 
is a solution if and only if (/?, 0, l^l) is a solution {\q\ := \J q^ + Q'^)- But we have 
seen that for arguments of the form (/?, 0, \q\), all the terms in equation ( [4.16| ) are 
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purely imaginary and hence we have reduced equation ( [4.16| ) to one equation in two 
unknowns. Accordingly define 

= (2cos/3 + (5i)| + -W{ip{0,p)+ij{(3,0,p)){l). 

We want to identify the set of zeros of the function g. Property (iii) of the last 
lemma implies that g{P, 0) = for all /?. We introduce the associated function 

9(P,P) ■■= • 

p 

Note that this function has a extension for p = 0. Furthermore 
ff(/3i,0) = ap5(/3i,0) 



(2cos/3i + di)l + -DWmdpi^ + = 0. 



We compute 



d/sg{pi,0) = dpdpgiPuO) 

= -sin/3i + -^DW{o)[dp^pfi]{i) + -^D^wio)[^p{^ + ^i),^/3^A{l) 

= -sin/3i/0. 

In the above calculation we used that DW{0) = and that dp^[l3i^Q) = 0, as 
//(/?, 0) = for all (5. Hence we can again apply the implicit function theorem and 
obtain in this way an even function I3{p), defind as map from a neighborhood of 
zero to a neighborhood of such that 

g{m,p) = ^- 

The evenness of /3 is a consequence of the uniqueness in the implicit function theorem 
and the evenness of g in p. We can summarize our considerations. 

Lemma 4.20 Given the function q) from Lemma 4-1% there exists an even 
function (3{q), mapping from a neighborhood of zero to a neighborhood of (3i such 
that the set of zeros of 

A{PMq)+QW{^{q)+^l{P,q))=0 
in a neighborhood Uq o/(/3i,0, 0) is given by 

q)eUo:q = 0}U q) e Uq : P = (5{^ql + ql)]. 
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Combining the last two lemmas we have proved the following theorem. 
Theorem 4.21 There exists a neighborhood Dq of in M^, a map 

r : Dq^ ii^u,, {qi,q2) ^ r{qi, 52) 

and a map 
such that 

(i) yqeDo: A(/3(Q))r(g) + W{r{q)) = 0. 



(a) r{q){—m) = r{q){m) for all m ^ "L^q ^ Dq. 

VgGA:||r((?)|k,, 
(iv) 3C > Oyq G Do : Eh>1 *Mk('7)MI < C\q\^- 

Proof : Let 

m ■■= m), 

r{q) := ^{q) + fJ-{l3{q),q). 

Then all the properties follow immediately from the last two lemmas. The differen- 
tiability of /3 at follows from the fact that /?(•) is an even function. 

O 



As we have proved Lemma 4.1£ and Lemma 4.20 using implicit function theorems, 



we can prove a uniqueness result for the constructed sequences r{q). 

Corollary 4.22 There exists a 5 > 0, such that the following holds. Let s G ^1^. 
and /3 E R, satisfying 

(I) < \\s\\i,u, < 6. 
(II) s(0) = 0. 



(III) s{—m) = s{m), for all m & Z. 

(IV) \p-/3i\<6. 



69 



y„,(t) := cn + ^ s(m) exp{im{f3n + 7^)), n Z, 

is a solution of 

yn = F{yn-i- Vn) - F{yn-yn+i), n E Z, 
then s = r[q) and (3 = f3{q), with 

q = {quq2) = 2{ Re{s{l)), /m(s(l))), 



and r, j3 are the functions defined in Theorem 1^.21 above. 



Proof : If (5 is chosen suitably small, it follows from Lemma [4.10| above that 



A{(3)s + Wis) = 0. 

Defining q as above, we see that the projection of s, Ps, satisfies assumptions (i) and 
(ii) of Lemma 4.1S| . Consequently we have Ps = n{j3,q). Furthermore the spatial 
frequency (3 satisfies the equation in Lemma 4.20| with g 7^ (otherwise // = s = 0, 



which contradicts assumption (I)). Lemma [4.20| yields that /3 = /3(g) which finally 
implies s = r{q) and the Corollary is proved. 



Remark 4.23 Heuristic explanation of the choice of the sign of (3. 

We investigate the transport of energy of the travelling wave solutions constructed 



in Theorem |4.21| . The amount of energy exchanged between the particles xq and xi 
is given by 



E 



F{xo{t) - xi{t))±o{t)dt. 



Note that E < means that on the average energy flows from xq to xi. Recall from 
equation ([4.2D in the introduction of this chapter, that 
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Let 



f{m) 
r'{m) 



(1 



)r{m). 



Then 



E 



27r 

7 



i'jmr{m). 



-r * r 



(0). 



Lfc=0 

Note from Lemma 4.19| (iv) that only r(l) and r(— 1) are of first order in q, all others 
are of higher order. We now evaluate the formula above. 

• 1 order: -^ao^'(O) = 0. 



2 order: 



27r 

7 



ai(f * r')(0) 



27r 

7 
47r 

7 



ai(f(l)r'(-l) +f(-l)r'(l)) 
ai Re [(l-e''^)|r(l)|2(-i7) 



47rai|r(l)p sin/3. 



• 3 order: no terms 7^ 0. 



In the system we are investigating we expect, that the driver excites outgoing waves, 
hence the energy should be transported in direction of increasing n. The above 
calculation shows up to third order in q, that this this is achieved for f3 with sin /? < 0. 
Therefore we have chosen solutions with /3 close to /3i rather than close to = — 



4.2 Construction of the periodic solution via Chap- 
ter 3 

The following theorem is the main result of the present chapter. 

Theorem 4.24 Let F, c, (6 m)rnGZj u! Satisfy the general assumptions and let F' (—c) < 
7^ < AF'{—c). Then there exists a neighborhood DofO such that for all e £ D there 



exist sequences n(e) G Cq^w;v{€) £ Ca^^w (o'o > defined in Definition 3.37 ) such 
that 



Xn{t) := cn + ebo - (u + w)(e)(0, 0) + ^ (u + t;)(e)(n, 771)6^"^^*, for n > 1, 

mez 

is a time periodic solution of the differential equation given by ( \3. j| j and (3.i). 
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Proof : We start the construction with r{q) E ii^w, which was obtained in 



Theorem 4.21 . Recall the definition in equation ( [4.5D , namely w{m) := w{m)(l + 
|m|)^, which again is an admissible weight function. According to the ansatz which 
was explained in the beginning of the present chapter, we define 



(4.25) 



u{q){n,m) := r{q){m)exp{iP{q)mn), for q G Dq. 



From Theorem 4.21 it follows that all of the conditions of Theorem |3.38 on the 
function u{q) are trivially satisfied, with the exception of condition (2), which states 
the differentiability with respect to the parameter g in a certain norm. It is obvious 
from Theorem 4.21 that each component of n is a function of q with 



_d_ 

dq. 



u{q){n, m) 



d d ~ \ ~ 

— — r{q){m) + ir{q){m)— — (5{q)mn ex.p{i(3{q)mn) 
oqj oqj I 



In order to show the dependence of u on g in the C-a^^w norm, we only have to 
verify that q -^u{q) is a continuous map from Dq into £„o-i,w It is straightfor- 
ward to check this by hand, using the following observations. 



exp{ip{q')mn) — exp(i/3(g)mn) 



< mn\ 



q\ sup 

sG[0,l] 



^PiQ + siq'-q)) 
dqj 



w{m)m?' < w{m). 
sup„>Q n^e~°"i" < oo. 



Therefore we are in a position to apply Theorem 3.38 of Chapter y and obtain 
this way the function v{q, e). Equipped with both, u and we turn now to Lemma 
3.31 of Chapter ^ Using Lemma 4.10| it only remains to solve 



(4.26) 



^(0, m) — n(0, m) = 0, for \m\ = 1, 



as by the conditions on 7 we have mo = 1- It suffices to solve equation ( 4.26 ) for 
m = 1, because for m = — 1 we obtain the complex conjugate of the equation. Hence 
there are two real equations to be solved (real part and imaginary part). Let 



g{q, e) := ehi - v{q, e)(0, 1) - u{q){Q, 1) G 
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Note that v{0, 1) is a function of {q, e) (see Theorem |]3| (v)) with \v{q, e)(0, 1)| < 
2(|e|+C|g|2) (see Theorem |3^ (ii), Lemma (iv)) and therefore ^^^(O, 0)(0, 1) = 
0. Furthermore ti(g)(0, 1) = ^{qi,q2)- We conclude 

9(0,0) = 

iDgg){0,0) = 

The imphcit function theorem then guarantees the existence of a neighborhood D 
of and a function q{e) such that g{q{e), e) = 0. This proves that for 




u{e) := u{q{e)), 
v{e) := v{q{e),e), 



the hypothesis of Lemma 3.31 is satisfied and we have successfully constructed the 



solution to the differential equation given by (B^) and (3.2). 



O 



73 



Chapter 5 

The Toda Lattice 



In this chapter we wih use the complete integrabihty of the doubly infinite Toda 
lattice {F{x) = e^) and show how the well known g -gap solutions contain a suffi- 



ciently large family of travelling waves to construct solutions of equations and 
(|3.2| ) for any number of resonances niQ. 

G-gap solutions were first constructed for the continuous analog of the Toda 
lattice, the KdV equation. Combining the complete integrabihty of the system 
as given in ([Q) with methods, developed in algebraic geometry (see 0]), these 
solutions can be expressed in terms of ratios of theta functions (see Q, p^M| ). 
Although these solutions have been studied and used in many different contexts 



(see e.g.[DT|, [ Krl |, [Kr2|,[McKT], etc.), we repeat their construction in the Toda 
case, because we will need some specific details about these solutions, which are 
not readily available in the literature. The resulting formulae will be evaluated in 
the case of small gaps and we derive a parameterization of time-periodic g'-gap 



solutions. This will be done in Section |5.1j . In order to keep the presentation self 
contained we provide the details of the construction in the Appendix^, where we will 
use only some general facts about hyperelliptic curves. The resulting parameters 
are closely related to those introduced first by T. Kappeler ([Kap|) for spatially 
periodic potentials in the case of KdV, then by T. Kappeler et al. [|BGGK| in the 
Toda case. It turns out that the condition of time periodicity determines the position 
of the midpoints of the gaps. In fact, there are two possible positions for each gap, 
corresponding to outgoing and incoming waves. We then proceed to prove in Section 
|5.2| that the basic result of Chapter ^ can be applied to obtain periodic solutions 
of the driven lattice for an arbitrary number of resonances. In the case of tuq < 1 



74 



these solutions are shown to coincide with those constructed for general lattices in 
the previous two chapters. Finally the essential spectrum of the corresponding Lax 
operator is investigated. Clearly it has a band - gap structure and we will determine 
the width of the gaps and their dependence on the Fourier coefficients of the driver 
to first order in e. 



5.1 The g-gap solutions 

First we briefly describe the construction of g-gap solutions via Baker- Akhiezer 



functions. We follow the construction in | Krl | (see also |]IM| ). 

Let g £ Nq = {n £ Z : n > 0} and denote by Rg the hyperelliptic curve of genus 
g, which is constructed by pasting together two copies of the Riemann sphere CU{oo} 
along the slits [£'0, Ei], [E2, E3], . . . , [E2g, ^29+1] ,where (£'0 < < • • • < £23+1)- 
Denote by vr the canonical projection of Rg on the Riemann sphere. We fix g points 
Pj ^ RgA<j < 9, satisfying TT{Pj) £ [E2j-i,E2j]. 

Then for each n S Z, t S R, there exists an unique Baker- Akhiezer function 
ipn{t, •) (up to multiplication by a constant) meromorphic on Rg\{Poo, P^}, with at 
most simple poles at Pj, 1 < j < g and a certain prescribed behavior at the essential 
singularities Poo and P^, depending on n and t. See Theorem |A. 11 in Appendix^, 
for the precise statement. The existence of these functions is proved by explicitly 
constructing them in terms of theta functions and uniqueness is a consequence of the 
Riemann- Roch theorem applied to Rg . Using the defining properties of the Baker- 
Akhiezer functions one is able to obtain functions Xn{t) such that the corresponding 
operators L and B (see ( |1.13| ) and ( 1.14 )) satisfy for all i G M and P £ Rg the 
following equations. 
(5.1) 



(5.2) 



(Li^Ut^P) = ^Tr{P)iJnit,P). 
^V'n(t,P) = (PV')n(t,P). 



From equations ( |5.1[ ), (5.2) we conclude for t G M, that 
(5.3) Lt = [B,L]. 

Hence Xn{t) is a solution of the doubly infinite Toda chain. Furthermore x„(t) can 
be expressed in terms of theta functions, namely 



(5.4) 



Xn{t) =nI + tR + fr (Un + Vt-Z) 
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The function fr : R^/Z^ — > R is essentially given by the logarithm of a ratio of 
theta functions with period matrix r. In general these solutions are quasiperiodic 
in n and t. The parameters /, R, r, U, V, Z depend on the spectral data chosen in 
the construction, i.e. they depend on the 3g + 2 real parameters Ei,0 < i < 2g + 1 
and Pj,l < j < g, or equivalently on {a,b, Xj,pj, Pj,l < j < g), where Xj denotes 
the midpoint and pj the half-width of the j-th gap. We are interested in obtaining 
the following choice for the parameters. 

(5.5) I = c. 

(5.6) R = 0. 

The next theorem will show that we can choose the parameters a, b, Xj,l < j < g as 
functions of the remaining 2g parameters pj, Pj,l < j < 5 such that equations 
dJ), dU), (13) are satisfied. 



Theorem 5.8 We assume that g G Nq, 7 > 0, c G R satisfy 

(5.9) 57 < 2e^t < (5 + 1)7. 

Then there exists a positive number 6 and functions a, b, Xj, Uj,T^^J^\ I < i,j < 
g, mapping {(pi, . . . iPg) G R^ : \pk\ < 5,1 < k < g} into R and which are even in 
each argument Pk,^ ^ k < g such that the following holds. 

For < pi, . . . ,pg < S the g -gap solution corresponding to the choice of param- 
eters 

Pj,Pj,a{pi, . . . ,Pg), b{pi, . . . ,Pg), Xj{pi, . . . ,Pg), 1 < J < fi-, 

is given by 

(5.10) Xnit) = era + m ^ , — — — 

^ ^ ^i-^U - Z\T)^{{n + ^)U + tV - Z\t) 

with 



(z) V 
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(ii) The map (Pi, . . . ,Pg) i — > Z is a surjective map from {S^y to R^/Z^ for all 
choices of parameters < pi, . . . ,pg < S. 

(Hi) 

(5.11) 7rir= diag [In pk) +T^''^^\ 



Furthermore formula \5.1C ) is also well defined for pj > and the function we 
obtain by letting some or all of the pj converge to agrees with the corresponding 
lower gap solution. 



Remark 5.12 (1) The proof of Theorem 5.8 is given in the Appendix ^ 



(2) In fact there are two possible choices for each function Xj The special choice 
made in the proof of Theorem |5.8| corresponds to the numerical observation, 
that gaps open up only in the lower half of the band. We will see in Remark 



5.30| at the end of the next section, that the physical reason for this lies in the 
direction in which energy is transported in the corresponding g -gap solution. 
The reader may recall that exactly the same situation occurred in Chapter 



with the choice of the spatial frequency /3 (see Remark 4.23) 



(3) The functions in Theorem 5.8 are also defined for negative values of the PkS. 
This extension is purely formal and is used to simplify regularity proofs at 
Pk = 0. 

5.2 Construction of the periodic solution via Chap- 
ter 3 

In this section we use the g'-gap solutions of the last section to construct periodic 
solutions of the driven lattice by means of the procedure in Chapter ^. Our goal is 
to prove the following theorem. 

Theorem 5.13 Let c,^, {bm)mez,w satisfy the general assumptions. Then there 
exists a neighborhood D of 0, such that for all e £ D, there exist sequences u{e) G 
^o,w,v{e) £ CfjQ^w (where ctq > was introduced in Definition 3.310 '^'^^ 



(5.14) Xn{t) := cn+ebQ-{u+v){e){<d,Q)+^{u+v){e){n,m)e:''^^\ forn>l, 



is a time periodic solution of the differential equation given by ( \3. j| j and (3.i). 
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Proof : By the general assumptions on 7, there exists a mo £ No, such that 

("^07)^ ^ ^ ^ {{mo + 1)7)^ 

We choose g := rriQ and thus satisfy the assumptions ( ^.9| ) of Theorem |5.8| . To put 
the results of the preceeding sections in a form that is suitable for the procedure 
of Chapter |3|, we still have to make some technical definitions and remarks. Let us 
begin with the definition of the parameters. Denote 

(5.15) pj := pj exp(27riZj) G C, 1 < j < g. 
The 2g real variables qj are then defined by 

(5.16) Pj = q2j-i + iq2j- 

Note that any choice of (7 in a sufficiently small neighborhood of in R^^ cor- 
responds to a choice of spectral data. In fact, let pj := \q2j-1 + iQ2j\- Furthermore 
Lemma B.47 in Appendix ^ shows that for any given choice of phase Zj G [0,1), 
there is a choice of points Pj which corresponds to the phase . Using ( B.63| ) we can 
now write equation ( ^.10 ) in the following form, 

(5.17) 40.W = .„ + ,„-i±I4?) +,„™i>. 

l + Tn+i{t,q) l + ro(0,g) 

with 

(5.18) r„(t,(?) := r{q){n,l)exp{-i{h+2l2 + ---+glg)-/t), 

«ez9\{o} 

(5.19) = J2 s{q){n,m)e'^"'\ where 




(5.20) exp (^27ri < I, U{p) > (n - ^)+ < /, r^'''^^) {p)l >^ 



(5.21) s{q){n,m) := Yl r(q)(n,l). 

ieZs\ {0} 

'1 H gig = -"1 

The Fourier series of the ^-gap solution is now given by 

(5.22) x^^\t) = cn+ ^ u(g)(n,m)e^'^^*, with 
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(5.23) u{q){n,m) 

(5.24) ui{q){n,m) 



ui{q){n,m) - ui{q){n + l,m) + U2{q)l{m=o}, 



-s{q){n,-THm) 



(5.25) U2{q) := ln(l + ^ r((?)(l, /)) - ln(l + ^ r{q){0,l)). 

iez3\{o} iez9\{o} 

The convergence of the above series follows from the smallness oi pi, . . . ,pg and will 

be verified below. 

Claim: 

There exists a neighborhood Dq of in M^^ and a constant C > 0, such that 



(I) for all q e Do : u{q) £ Co,w and ||'u(g)||o,«) < C\q\. 
(II) for all qGDo: E|m|>n^o ^Ml^('7)(0, HI < C\q\^- 

(III) u : I?o ~^ C,-(Ti,w, Q ^ u{q), is a map, where ai was introduced in Definition 
[3:37 . 

(IV) f : Dq ^ R^^,g 1-^ (n((7)(0, m))^^]^ is continuously differentiable and 

det(i?,/(0)) / 0. 

(V) xi°^(t),n G Z is a smooth, real valued solution of the doubly infinite Toda 
lattice. 

Before we proceed to check all these properties, let us show that they suffice to 
prove the theorem. First we have to verify that u{q) satisfies the conditions (l)-(3) 
in Theorem p.38| , but this is an immediate consequence of (I), (III) and (V). Hence 
we obtain a v{q,e) G C^q^w, satisfying conditions (i) to (v) of the Theorem 3.38| . 
Using (II) we obtain in addition, that 



(5.26) 



\v{q,e) 



<2{\e\+C\q\^). 



We now show that it is possible to satisfy (1), (2) and (3) in Lemma 3.31 

(1) follows for u{q) from (V) and the Remark |3.32| . 

(2) follows for v{q, e) from Theorem ^.38| . 

By Theorem p. 38 it suffices to solve (3) for 1 < m < g. This gives 2g real equations, 
for which we have the 2g real variables qj available. Define 



9{q, e) := i^bm - u{q){0, m) - v{q, e)(0, m)) 



m=l 



?29 
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Using Theorem 3.38| (v), equation ( ^.261) and what we know about function / from 



claim (IV), we see that g is a function and 

5(0,0) = 0. 
D,g{0,0) = -DJ{0), 

which is an invertible matrix by (IV) . The imphcit function theorem yields a function 
q{e), such that g{q{e),e) = 0. Thus 

u{e) = u{q{€)). 
v{e) = v{q{e),e). 

satisfies all the conditions of Lemma |3.3l| . 
It remains to prove properties (I)-(V). 

(I): 

It is easy to read off the following estimates from the above definitions. Denote for 



3C,(5>0:V|g| < 6,n e No,l e : \r{q){n,l)\ < Cl'l'l^ll^l". 
/i + 2/2 + • • • + gig = —m implies that |/| > ■^3-. Furthermore 

3C,d>0:y\q\<6,neNo,kQeN: ^ |r(g)(n, /)| < Cl''ol'|g|l'=ol'. 

\l\>ko 

Together with the observation, that the sum defining s{q) does not contain the term 
where Z = 0, it follows, that 

2 

3C,<5 > : V|g| < 6,n £ No,m £ Z : \s{q){n,m)\ < (C|g|)"'^''^^'^\ 
Because of the estimate on the weightfunction w{'m) < Ce'^™', for some constants 



C,a > (see Definition 3.21 and below), we conclude that s{q) G £0,10 for \q\ small 
enough and 

3C,6>0:y\q\<d: \\s{q)\\o,^ < C\q\. 

Equation ( ^.24]) shows, that ui{q) G £0,10 and ||Mi(g)||o,«; < C\q\ for \q\ small enough. 
Applying the same kind of estimates again we obtain 

3C,5>0:y\q\ <5: \u2{q)\ <C\q\. 
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This completes the proof of (I). 



(II): 

The definitions at the beginning of the proof yield 

k 



V|m| > mo : n(0,m) = s(0, m) - m) + ^ — j- (^s(0, - s{l,-Y^{m)j . 



k=2 



Reworking the corresponding estimates in the proof of (I) one obtains 

2 

3C,(^>0:V|g| < 5, nE No, |m| > mo : \s{q){n,m)\ < (C|g|)'"^''^^'^\ 
from which we easily conclude that 

3C,(5>0:V|g| <5: ||s(^7)l|H>„„}||o,»<C|g|2. 
On the other hand the estimate above on 115(9)110,10 implies, 



3C,6>0:y\q\<6,neNo : II E 5(g)(n, •)*'=||o,«, < CjoP 



k=2 



This proves (II). 
(Ill): 

Let us start with the differentiability of each r{q)(n,l) with respect to Qj. From 
equation ( 5.20| ) and Theorem f).8\ we learn that the only possible problem lies at 



points where one of the pj = 0, as pj = y Q2j_i + l2j ^'^^ ^ differentiable function 
of q in those points. However, the functions U and t^^^^^ are even in each variable 



Pj (see Theorem 5^, lj{lj — 1) is even and pj = (?2j-i +iQ2j- Hence r{q){n,l) is C^. 



Simple estimates show that 

3C,C,5>0:y\q\ < 5,n£ No, I eZ<^\{0} : \—r{q){n,l)\ < C'Cl'l'|g|l'l'~^(n + 1). 

oqj 

3C,C,6>0: V|g| , \q'\ < <5, n G No, / G \ {0} : 

\A.riq')in,l) - -^r{q){n,l)\ < CCl'l' max(|(?'|, |g|)--(l'l'-2'0)(n + l)2A(g', g), 
with 

A(,'.,) := max L' - ,|. | - ||-.("«)(,') - |-r(-)(,)|) . 
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Note that the powers of \l\, which are produced by the differentiation have been 
subsumed into the C'''^ term, simply by increasing the constant C. 

Our next goal is to prove that q ^ s{q) is a map into C-a-2,w, with a2 '■= 
The differentiability of r{q)(n,l) with respect to q and the absolute and uniform 
convergence of the sum 

Z G Zf \ {0} ^ 
h -\ gig = -m 

shows that for all n G No,m G Z : s{q){n,m) is a function of q. Therefore we 
have only to prove that the map q ^ -£ps{q){n, m) is a continuous map into C-ai^w 
This, however, follows by applying the arguments given in the proof of (I) to the 
estimates for and |^r(g')(n,/) — ^r(g)(n,/)| given above and from 

the observations that A(q', q) ^ 0, as q' ^ q and that sup„>o('^ + l)^e~'^^" < oo. 

Let us now investigate the differentiability of ui{q). Employing the calculations 
we made in the proof of Theorem 3.38| , step 4, where we were in a similar situation. 



we obtain 

^m{q) = '£{-l)'^s{q){n, * T^s{q). 



Standard arguments yield the continuity of the map q ^^i(^) ^~2a2,w = 
r 

The differentiability of the function U2{q) can be established from the differen- 
tiability of r{q){n,l), the absolute and uniform convergence of the corresponding 
sums (see equation (|5.25| )) and the differentiability of the logarithm away from 0. 
The proof of property (III) is completed. 

(IV): 

The differentiability was already proven in (III). We only have to investigate the 
terms of first order in q of u{q)(0,m) for 1 < m < g. Following the reasoning 
in the proof of (II), we see that we only have to consider the first order terms of 
s{q){0,m) — s(g)(l,m) which are given by r{q){0,em) — '''{q){l,em), where Cm '■ = 
(0, . . . , — 1, . . . , 0) denotes the unit vector in with —1 at the m -th entry. This 
implies 

(5.27) VI < m < 5 : u(g)(0, m) = p„2isin(7r[/™(0)) exp(r4^f^)(0)) + Odgj"). 
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We are done if we can show that sin(7r[7m(0)) 7^ 0. It turns out that we can compute 
this quantity exphcitly. By Lemmas |B.19| , |B.36 and Proposition B.62 we conclude, 
that 

(5.28) Um{0) 
Therefore 



2 ^ b-\m 
- — arctan \ / 



sin(7rC/^(0)) 



-2 tan arctan 



' b- Xm 
Am. — a 



(5.29) 



b — a 
2 

b — a 



\J {Xm - a){b - Xm) 



/ 0, 



by Lemma B.41 



cos arctan 



Xm — , 



(V): See Appendix 0, Theorem [AI21 



O 



Remark 5.30 The choice of Xj (compare with Remark 5.1^ (2), see also ( \B.(\ )). 



In Appendix ^ it is shown that the condition of time periodicity of the g-gap solution 



leads to equation (B.42) in Lemma B.4L This equation was solved by choosing the 
Xj 's as functions of the remaining parameters. As it was stated in Remark ^.12| (2), 
there exist two choices for each Xj, one in the lower half of the spectrum and one 
in the upper half of the spectrum. We will now discuss the difference of these two 



choices. Proceeding as in Remark 4.23, we investigate in which direction the energy 



is transported in the g-gap solution up to second order in q. The arguments given 
in the proof of property (IV) in the proof of Theorem 5.13| above yield 



X 



(0) 



it) 



me 



iymt 



+ 0{\q\ 



m=l 



c+ E n((7)(0,m)exp(-27ri[/„(0))e'^'"* + O(|<7| 



m=l 



where u{q){0,m) was determined in ( 5.27| ). Now we compute the energy which 



is exchanged between particles x'n^ and x^'' during one period. Repeating the 



(0) 
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calculations in Remark [4.23 we arrive at 



E = - r F{x'^>{t)-x';'>{t))x'^>it)dt 



9 



47rexp(-c) ^ m|u(g)(0, m)!^ sin(-27rC/„(0)) + 0(1^1 



m=l 



Each term in the sum corresponds to the energy transported by one phase of the 
multiphase solution. It is transported in the direction of increasing n, if 

sin(-27rC/„(0)) < 0. 



Equation ( p. 28 ) implies that 



Xrr 



-2iTUmiO) = 4arctan 



Therefore -27rC/^(0) G (7r,27r), if Ar„ < ^ and -27rC/„(0) £ (0,7r), if A„ > 
Thus the choice we make in ( p3.6| ), corresponds to a solution where the energy is 
transported outwards in the direction of increasing n. As described above our choice 
implies that all gaps open up only in the lower half of the spectrum. 

Remark 5.31 Comparison with solutions of general lattices for g = 0,1. 

Recall, that for mo = and mo = 1, we were able to construct periodic solutions 



for general lattices (see Sections 'iA and ^]2| respectively). We will verify that these 
solutions are the same as we have constructed in the present section in the Toda case, 
by showing that they produce the same families of sequences of Fourier coefficients 
u{q), which are used in the construction described in Chapter ^. 

The case g' = is trivial, as the 0-gap solution is simply given by = cn, i.e. 
u{q) = 0, which is the choice we made in Section 3.4. 



In the case g = 1 equation (5.10) and the periodicity of the theta function in the 
real direction shows that the one-gap solution is of a form to which Corollary 4.22| 
applies. We have to check the assumptions (I)- (IV) in the Corollary. (I) and (III) 
are trivially satisfied. 

Property (IV) is satisfied, if we can prove, that —2ttU{0) = /3i, where /?i was 



defined in ( |3.9| ). By the Remark 5.3C above we know already that both quantities lie 
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in (vr, 27r) (mod 27r) and therefore it suffices to show that cos(— 27r[7(0)) = cos(/3i). 
Equation (^) yields cos(/3i) = — Using ( 5.29| ) we compute 



cos(-27r;7(0)) + ^ 



872 



(6 - ay 



+ 



-1 + 



7 



2 exp(— c) 



as by the choice of b^'^^ (see ( B.55| )) in the proof of Lemma B.52, we have b — a 
4exp(-§). 



Property (II) of Corollary 4.22 is not satisfied as the zero Fourier coefficient 

s(0) = In. 



^{^U-Z\t) 



^-^U-Z\t) 

does not necessarily vanish. But this corresponds only to adding a constant to 



u{q){n,0), which does not change the equation for v in Lemma 3.31| (note that 



W{u, v) only contains Au, see Definition 3.2§| ) and which does not change the value 



of a{n,m) in Lemma 3.31. 



Remark 5.32 The opening of the gaps to first order in e. 

We will now investigate the size of the instability regions and their dependence 
on e and the Fourier coefficients of the driver. This question of basic interest was 
well studied in the continuous case for periodic potentials (see |MW] and references 
therein). In our situation recall that L denotes the semi- infinite Lax operator corre- 
sponding to the solution constructed in Theorem 5.13 above (see also 1.27| ) . As L is 
obtained from the doubly infinite operator L by restricting it to n > 1 and adding 



an operator decaying exponentially in n (which corresponds to the v-term in 5.14) 
standard arguments of spectral theory imply that 

o"ess(-^^) = (^essiL), 



and hence 

(See equation ( |5.1D ). Therefore the width of the j-gap in the spectrum of L{t) is 
given by \pj\, 1 < j < g- They were determined as functions of e, when we solved the 
resonance equations ebm = u{q){0,m) + v{q,e){0,m), 1 < m < g in Theorem 5.13| . 



'a Xi — Pi 


U 


'Xi+Pi \2-P2 


U • 


• u 


'Xg+Pg 


6" 


2' 2 




2 ' 2 . 






[ 2 


2 
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As v{q, e) is a fixed-point of the operator T{q, e, •) given by ( 3.45| ) in Chapter ^, we 
conclude from Theorem |3.38| (ii), Proposition |3.29| (ii) and the fact that q{e) = 0(e) 
that f (g, e)(0, m) = O(e^), for 1 < m < g. Equation ( p.27| ) then yields that the 
following relation holds. 



(5.33) 



|e6^| = 2| sin(7rf/^(0))| exp{T^2\0))\Pm\ + 0{e^). 



Below we will show that 
(5.34) 



T-('-e9)(0) 



In 



8(m7)^ 
b — a 



Using (|5.2g| ) and ( |5.33| ) we can express \pm\- 

(5.35) \pm\ = 2\ebjn\m'y + O(e^) 



We compare this formula with numerical results shown in Figures |l.3l| and |l.32| . We 
recall the choice of parameters, \ebi\ = 0.1, |e62| = 0.05. Up to first order, equation 



( 5.35 ) yields 

• 7 = 1.8 :pi 0.36. 

• -f = l.l:pi^ 0.22, p2 ~ 0.22. 

These values are in good agreement with the numerical experiments. In fact. 



7 = 1.8, Figure 1.31 : pi = 0.34 ± 0.007. 



• 7 = 1.1, Figure 1.32 : pi = 0.223 ± 0.005,p2 = 0.215 ± 0.15. 



We will now sketch the derivation of formula ( ^.34 ). By Theorem 5.5 and Proposition 



B.62| (vi), we may consider the one-gap situation, where only the m -th gap is open 



and all other gaps are closed. Using equations (1R32|) , (p^ , (p^ (^1211 ), (pl^) 



and (B^) and Lemmas B.14 , B.IS we arrive at 



T-{»-e9)('0) 



- — lim 



— ,, -.dE — hm In — 

^{E-a){b-E)[{E-X^Y-pi,] ™ 



(5.36) 

The quantity hm was defined in ( [B.13| ) and is given in this case by 

1 



(5.37) 



^J{\ra - a){b - Xra) 
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We define the auxiliary functions 

(5.38) f{E) : 

(5.39) giE) : 
Then 



^iE-a)ib-E) 
f{E) - fjXm) 

E — Am, 



Cr^rn—Pm 1 \ 

/ — ,, zdE — hm In — 

Ja ^{E - a)ib - E)[{E - - pI] Pm^ 

lim I r"^'""^ f{^m) ,^ _ /^^ In . ^ ^ 



Pm->0 y/{E - Xraf -pin 

+ lim ( r-"- 4E=M$B=dE] 

(5.40) = (/) + (//). 

Using the appropriate changes of variables we evaluate 

(5.41) (/) = /i„ln2(A^-a). 

At, 



{11) = / g{E)dE 

J a 



(5.42) = hm I / ' ^^^] dE - In -^"^ 



The remaining integral can be integrated and we obtain 

(5.43) r-^-^ 4^dE = ^ f m i + m + ln(l + 0(.)) V 
^ ' Ja ^ - Am e b-a ^ ^ 

Note that for Am = Xm{Pm = 0) equation ( [B.43D yields hm = Therefore we can 
determine 

(5.44) (//) = /i„fln4^-ln(Am-a)V 



(5.45) (I) + (11) = hmln 



b — a 

8(771,7)^ 
b — a 



It follows from (Q) and (|^ that r^'fm (0) = -In^^^lL^ which proves ( jO^ ) 
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Appendix A 

Definition of the g-gap solution 



It is our goal to derive a formula for the well known (?-gap solutions of the To da 
lattice. 

Let g £ No = {n £ Z : n > 0} and denote by Rg the hyperelliptic curve of genus 
g, which is constructed by pasting together two copies of the Riemann sphere CU{oo} 
along the slits [Eo,Ei], [^2, -^3], . . . , [£^29, -^29+1] ,where (£'0 < -Ei < • • • < £29+1)- 

A point on Rg is denoted by P, and the canonical projection of Rg on the 
Riemann sphere is given by vr. We write E = 7r(P). 

• For 1 < k < g let Pj £ Rg be a point in the j-th gap, i.e Tr{Pj) G [£'2^-1 , E2j]. 



ak,l3k, ^ ^ k < g denote the canonical homology basis for Rg (see [FK, III.l] 
for a definition of a canonical homology basis) as shown in the Figure |A.l| 
below. 




Figure A.l: The canonical homology basis of the Riemann surface 



The existence and uniqueness of the following differentials can be deduced from the 



Riemann Roch Theorem (see ||FK| , III. 4. 8]) applied to Rg (see e.g. pS| , VI. 4, Satz 
32], Remark 16.17]). 
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Lo := {i-^i)i<i<g is defined as the unique basis of holomorphic differentials on 
Rg, whicli is normalized by the condition /^^ loi = d^^i- Note that all the 
holomorphic differentials on Rg can be written in the following form. 

2g+l 

(A.3) RiE) := Yl{E-E,), 

3=0 



where p is any polynomial in E of degree < g — 1 (see [FK, III. 7. 5, Corollary 
!])• 

is the unique meromorphic differential on Rg with simple poles at Poo and 
(residues: 1 /resp. —1), holomorphic everywhere else and normalized by 
/^^ uj^^'^ = 0, for 1 < k < g. Here the points Poo and denote the point at 
infinity on the upper and lower sheet respectively. We can write 

where p^^^ denotes an uniquely determined polynomial of degree g. 

the unique meromorphic differential on Rg with second order poles 
at Poo and (principal parts: i^d^ /resp. ^d£,, where ^ is the coordi- 
nate around oo : ^ := :;;^py), holomorphic everywhere else and normalized by 
u;(2) = 0, for 1 < < g. Again this differential can be written in the 
following way. 

where p^'^^ denotes an uniquely determined polynomial of degree g + 1. 
We define the following /3 periods. 

(A.6) T := iTk,i)i<k,i<g, with r^,/ := / w;. 

The Riemann bilinear relations imply that r is a symmetric matrix with pos- 



itive definite imaginary part (see |FK, III.3.1, III.3.2]). 
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• We denote K to be the vector of Riemann constants with respect to the base- 
point E2g+i (compare |F^, VI.2.4]). 

The following three integral functions are multivalued and depend on the path 
of integration. 

MP) :=/4..- 

0(2)(P) :=/4^^^(2) oni?,\{Poo,Pi}. 

We recall the definition of the Riemann theta function (see e.g. [FK, VI. 1]). For 
a symmetric matrix with positive definite imaginary part, 

where < u,v >:= J2j=i'^j^j- The "periodicity" properties of the Riemann theta 
function are as follows. Let 1 < m < g and denote by Cm the m -th column of the 



g X g identity matrix and by Xm the m-th column of x- Then (see e.g. [FK, VI. 1.2]) 
(A.9) Hv + em\x) = ^vlx)- 

(A.IO) '&iv + Xm\x) = exp{-27rivm - mXm,m)'&{v\x)- 

The main tool for the construction of the g-gap solutions is the following exis- 
tence and uniqueness theorem for the Baker - Akhiezer function. Using the above 
notation, we have: 

Theorem A.ll (^r^, see also /p^j For all n G Z,t £ R, there is a 

unique (up to multiplication by a constant ) function ipn{t, which is not identically 
equal to and satisfies 

(i) tpnit, •) is meromorphic on Rg \ {Pao^P^}- 

(a) ipnit, ■) h<is at most simple poles at Pi, . . . ,Pg £ Rg\ {Poo,P^} o^nd is holo- 
morphic elsewhere on Rg \ {Poo, P^}- 

(Hi) tpnit, ■)E^"'^xp{±^E) has a holomorphic continuation at Poo /resp. P^. 

Proof : This result is well known and we only sketch the proof. First we prove 
existence. Define 

(A.12) *„(t,F) := e.p(„n<')(P) + '"'^'(P)) '"-'"',',^,^; ! ^'.^ 
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where n^^\n^^\A, U, V and r were defined above and Z will be introduced below. 
Before we can check that this function satisfies all the conditions of the theorem, 
some more remarks are needed. 

• As already mentioned above the functions A,VL^^\Q.^'^'^ are multivalued. We 
will show below that nevertheless the function is well defined if we 
insist that the path of integration is the same for all three functions. 

• Z := K + ^i<m<g -^{Prn) ■ As we Want ^n(i,-P) to be well defined we now 
specify the path of integration from -E'2g+i to P^- We first integrate on the 
upper sheet along R on the C+ := {E G C : Im(E') > 0} side from £"29+1 to 
the branchpoint E2m-i and then from £^2m-i to Pm along the real axis, where 
the path always first stays on the upper sheet, and in case that lies on the 
lower sheet, we switch the sheet at E2m- Of course we could have chosen any 
fixed path from -^29+1 to Pm-, but the preceeding choice leads to an especially 
simple formula for Z in ( B.46| ) below. 



• In the case that 5 = the theta functions are simply replaced by the factor 1. 
One checks that all the steps of the proof given below are trivially satisfied in 
this case. 

In order to see that ^n{t^ P) is well defined we have to examine what happens if 
we add to the path of integration one of the cycles a^, 1 < A; < (7 or 7^, j = 1, 2, 
where the 7 's are the cycles around Poo /resp. P^. The 7 cycles only effect 
$7^^) by adding a multiple of 27ri, which does not change the value of ^'„(t,P) 
because of the exponentiation. The a cycles only change the entries of A{P) by 
adding integers which has no effect because of the periodicity of the theta function 
in the real direction (see ( A.9D ). The /3 cycles change all three integrals and it is 



straightforward to show that all the factors cancel out by the monodromy property 



of the theta function ( A. 10 ) and by the choice of the vectors U and V (see ( A. 7 ) 
and (ra). 

are holomorphic in Rg \ {Poo;-Pi>}, properties (i) and (ii) are 
equivalent to showing that the zeros of '&{A{P) — Z\t) are all simple and that they 



are given by Pi, P2, . . . , Pg. By definition of Z and by |FK, Theorem b,VI.3.3], 
(see also [|G1| , Thm 17.9]) this is equivalent to proving that the divisor P1P2 ■ ■ - Pg 
is nonspecial, i.e that there exists no holomorphic differential on Rg vanishing at 



91 



Pi, P2, . . . , Pq and which is not identicahy equal to zero. But this is easily deduced 



from the characterisation of holomorphic differentials as given in ( A.2 ) and from the 
fact that vr(Pj) 7^ T^{Pj) for i 7^ j. 

Property (iii) can be verified from the definitions of 0;^^-* and uj^'^\ 
This settles existence and we can now turn to the question of uniqueness. Sup- 
pose ^n{t,-) satisfies (i),(ii) and (iii). We consider the function . It is mero- 
morphic and its poles are the zeros of '&{A{P) + Un + Vt — Z\t). It will be shown 
below (see Remark p3.48 ) that for all n and t, 'd{A{P) + Un + Vt — Z\t) has exactly 



g zeros, one in each gap, which form a nonspecial divisor and by the Riemann-Roch 

in(ir) 



theorem it follows that ?"[|''| must be a constant (see e.g. [ [FK] , III. 4. 8], |G1, Thm 



16.11]). 

O 

We now introduce a normalization of the Baker -Akhiezer function. Denote by 
ipn{t,-) the uniquely defined B A- function which has the following Gxpansion at Poo- 

t ^ 

(A.13) i;n{t,P) = E-^exp{--E){l + Y,it{n,t)E-'). 



This is possible as we see from ( A.12| ) and the position of the zeros of i?(j4(-) + 



Un + Vt — Z\t) described above, that ^n(i, -foo) 7^ 0. The expansion at P^ is then 
written as 

t °° 
(A.14) Mt,P) = E^expi-E){(^in,t) + Y.Cin,t)E-'). 

We can express ^'^{n,i) in the following way. Let us expand ^'„(t,P) as given 
in equation ( A.12|) around P^ and P^. From the definition of lo^^^ and cj*-^-* it is 



obvious that we can expand their integrals in the following way around the infinities. 

00 

(A.15) ^^^\P) = ^(In^; + ^I;^^"') around Poo /resp. P^. 

z=o 

00 

(A.16) rj(2)(P) = ^-{E^Y^ RfE~^) around Poo /resp. P^. 

^ /=o 

In order to have the zero order terms Iq and well defined we shall now fix 
the path of integration from E2g+i to Pqo and P^ as the path along the real axis 
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from E2g+i to +00 on the upper sheet /resp. on the lower sheet. The holomorphic 
continuation of "^n{t, P)E^Gxp{^) takes at Pqo the value 



2 fl(,4(P„) - Z|t) 

The path in the integration of to is chosen as above, which determines the value 



of A{Poo) uniquely. It follows by the Riemann bilinear relations (see e.g. [FK 
111(3.6.3)]), that 

(A.17) ^ f = r^u; = A{P^) - ^(P4). 

2tti Jp Jp^ 

As u> differs on the different sheets only by a sign we obtain with ( [A.7| ), that 
(A.18) A{Po.) = ^U = -A{P:,). 

At P^ the holomorphic continuation of ^n{t,P)E~'^exp{—^) takes the value 



>)(/l(Pi) - Z|t) 



Hence 



(A.19) 



t 



•-0 



^{U{n - i) + Ft - Z\t)^{W - Z\t) 



d{U{n + \) + Vt- Z\T)d{-\U - Z\t) ■ 

In Appendix ^ we will determine formulae for Iq , Iq ,Rq , Rq , U,V,Z, from which 
it is easy to see that all the parameters are real and hence ^(^(n, t) G M \ {0}. But 
(0, 0) = 1 and by the hence by continuity (regard n as an arbitrary real variable) 
^o{n,t) > 0. We define 

(A.20) Xn{t) := In ^^{n,t) 

as a real number. The next theorem shows that Xn{t) is a solution of the Toda 
lattice. 



Theorem A.21 (^Kr^, jgl/ j Let x„(t) be defined as in l[A.2(\ ). Then twice 



differentiahle (in fact infinite differentiahle) for all n £ Z and furthermore 

VnEZ,tGM : Xn{t) = exp{xn-iit) - Xn{t)) - exp{xnit) - Xn+l{t)). 
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Proof : Recall from the Introduction that we can write the Toda equation in 
Lax pair form by using Flaschka variables (see ( |1.12 ) - ( 1.15D ). We will show below 
that for 
(A.22) 

the following equations are satisfied for all t G M and P & Rg\ {Poo, P^}- 



(A.23) 



(A.24) 



Li, = ^Ei,. 
1^ = ^^. 



Equation ( |1.15| ) follows as the compatibility condition of equations ( A.23| ) and 



(|A^). 

Proof of equations (|A.23D and ( [A.24| ). 

We define the auxiliary function Ai := {{L — ^)ip)n- Using ( A.13| ) we can calculate 
the behavior of AiE"-exp{^E) around Poo and we arrive at 

^exp{-^){Ctin - 1, i) - Ctin, t) - Xn) + O(i^-i). 



In the same way we determine from equation ( A. 14 ) the behavior of AiE ^exp{—^E) 
around P^. It is given by 

ie^p(^)(e-n+i^-(^ + 1, t) - e-^"er(n, t) - x„) + 0{E-'). 



The quantity Ai satisfies conditions (i) and (ii) of Theorem A. 11 and thus we can 
conclude that 



(A.25) Clin - 1, t) - t) - x^ = e-^"+'Ci{n + l,t)- e-""er(n, t) - x^. 



Proceeding in an analog way for the auxiliary function A2 = ((^ — B)i))n we 



conclude again from Theorem [A.ll| that 



(A.26) Ctin-ht)-Ctin,t) 



-a;n+i 



Adding equations ( A.25|) and ( A.26|) we see that Cf{n — l,t) — C,i{n, t) — Xn = 
and again by Theorem A.ll| it follows that Ai and A2 are both equal to zero and 
this proves equations ( A.23| ) and ( [A.24 ). 



O 
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Appendix B 



Evaluation of the g-gap 
solution for small gaps 

The goal of this appendix is to prove Theorem |5.8| , that is to produce a 2g-parameter 
family of time periodic solutions of the doubly infinite Toda lattice with period ^ 
satisfying certain regularity conditions. We use the g -gap solutions produced in the 
last section in the case that the gaps are small. First we introduce some notation. 
We choose the variables Xj,pj,l < j < <? in such a way, that -£"2^-1 = — Pj 
and E2j = Xj + pj,l < j < g. In addition we call a := Eq and b := £"23+1. 
Hence the g-gap solutions are completely parameterized by the 3(7 + 2 real quantities 
a,b, Xj,pj, Pj,l < j < g. The plan of the present Appendix is as follows. We will 
determine the dependence of the quantities t,U,V, Z, , in equation ( [A. 19 ) 



on the parameters a,b,Xj,pj,Pj,l < j < g. Then we will proceed to show that 
for small gaps there is indeed a choice of parameters, such that the corresponding 
solution Xn{t) as given by ( A.20| ) has all the properties of Theorem |5.8| . Special 



emphasis will be given to the limits pj — > 0, i.e. when gaps close. Even though the 
analytic expressions will achieve some limit value it is not clear a priori, that this 
limit coincides with the formula for the corresponding lower-gap solution. Therefore 
this has to be checked separately. 

Remark B.l During the Appendix we will always assume that the hypothesis of 



Theorem 5^ holds, i.e. 

57 < 2e~t < (5+1)7- 
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We specify the range for the different parameters. 



(B.2) a G (-2e-f -ea,-2e-f +ea). 

(B.3) b G (2e-t -ef„2e-t +eb). 

(B.4) p, G (0,ep^). 

(B.5) Xj G ( Af - eA, , Af + eA, ) , where 

(B.6) Af := ^-^(^) -.V,l<.<.. 

(B.7) ^(Pj) G [Xj-pj,Xj+pj]. 

The main requirements which have to be satisfied for the choice of the various e 's 
is that we have to ensure that the bands do not vanish, i.e. a < Ai — pi < • • • < 
Xg + Pg < b. Furthermore we want Af to be real, which can be achieved by the 



assumption (5.9) above together with the choice of a and b. During the calculations 
some other conditions on the smallness of the e 's will occur (e.g. induced by the 
use of the implicit function theorem) and of course we want them to be satisfied as 
well. 



B.l The holomorphic differentials and the r matrix 



The holomorphic differentials on Rg can be written as -^i^=dE (compare with 



A.2), where p is a polynomial of degree < 5 — 1 , 
29+1 



R{E) =\{{E- E,) = {E-a){E- b) \{{{E - X,f - p] 



^1) 

j=0 j=l 



and y/R{E) is defined on Rg with the usual convention, that on the upper sheet 
y/R{E) — > +00, as E ^ +00. Our first goal is to determine the canonical basis 
of the holomorphic differentials uJi = --jiS]=dE, which has to be chosen such that 

la, = 



k,l- 



For the calculations it turns out that the following basis of polynomials of degree 
< g — 1 is useful. We denote 

(B.8) e,{E):= f[ ^ " 



Xj — Xr 
m = 1 
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We begin with the calculation of some elementary integrals. The first candidate 

is 

and 

(B.IO) B := {h,j)lj^,. 

Here the square root in the denominator always denotes the positive root. In order 
to demonstrate how this integral will be analyzed (especially in the limit as pk — >■ 
0) we split the integral in two parts. Let d be any fixed number in the interval 
[E2k-2, E2k-i], which is away from the edges of the bands independently of the 
choice of the parameters as they vary over the allowed regions. Let 

m = 1 
m ^ k 

then there exists a smooth function fj^k such that we can write 
(B.12) fj,k{E) = 4fe(Afe) + {E- Xk)fj,k{E) 

for E G [d, Ajfc] . By simple changes of variables we have 



1 ds 

zdE 



'Afc - d 



= arccosh 



Pk 



, ..Afc - d , / 1 /l pi 
\n2— hln - + ' ^ 



Pk 12 V 4 4(Afe-d)2 



Note that the last quantity is a function of all the parameters, namely (a, b, \j,pj). 
Furthermore we see that the expression just depends on pj and therefore we can 
define it for pj < such that the smooth dependence on the parameters is preserved 



97 



and we obtain a function, which is even in all (small) pj 's. This continuation is 
purely formal. Using the abbreviation 

(B.13) hk := TT 



V(Afc - a)(b - Afc) y^{Xk - Xm)'^ -pV 

m = 1 
m k 

we state 

Lemma B.14 Given the definitions as above, then 

(B.15) bkj = Skjhk In — + 4_i In — h Wf\ 

Pk Pk-i 

where b^^Y^ ^-^ ^ function of the parameters {a,b, Xj,pj), which can be extended 
for nonpositive values of the pj 's in an even way, preserving the smoothness. 

Proof : We split the integrals as described above. The integrals from E2k-2 
to d can be dealt with in a similar way. 

O 



There is a second elementary integral which we wish to discuss. Define 

Utilizing the notation introduced in equation (|B.11[) we obtain without effort 

/T^1'7^ ~ fjMM+Pks). 

(B.17) ak,=j_^ ds, 

from which it is obvious that dkj is a smooth function of the parameters and can be 
extended to (small) nonpositive values of pj as an even and smooth function. From 
the formula ( B.IT] ) we can also read off that 



(B.18) Ofcj = (5fcj7r/ifc + a-^"^ 



with d^i^j^ = 0{p\).{ (ho) stands for higher order). Let us introduce the matrices 
^ := (^fcjOfc j=i) A^^°^ := (a[,^°^)^ j-^]^ and C := . We can again state a lemma. 
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Lemma B.19 Given the definitions as above. Then 
(B.20) A = diag (nhk) + A^''°l 

(B.21) 



C= diag (J_) + 



All entries of the matrices are smooth in the parameters {a,b, Xj ,pj) and have an 
even and smooth extension for (small) nonpositive pj . Finally we have the estimates 



yi<j,k<g: aff = 0{pi) and c^J = 0{pi). 

Proof : The claim has been shown for A. Expanding C in a Neumann series 
completes the proof. 

O 

Using these two lemmas we can evaluate the corresponding integrals over the 
cycles of the canonical homology basis. Let 



(B.22) 



(B.23) 



^^^dE and A:=(afc,,)f . 



^^^.dE and B := {bk,j)i 



Looking at Figure A.l in the beginning of Appendix |^ and keeping track of the 



signs we obtain the following formulae 

(B.24) A = diag {2{-l)3-^)A. 



I 1 



(B.25) 
with 
(B.26) 
(B.27) 



B 



diag (2i(-l)s~'^)B. 



V 1 ■■ 1 



diag {-2i{-iy-''hklnpk), 
/ 1 \ 

: •., diag {2i{-l)s~^)&"'^\ 

[ i ... I J 
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We return to determining the basis of the holomorphic differentials 

n{E) 



(B.28) 



-.dE. 



Vim 

Obviously we can write ri{E) = X]f=i '^i{^j)^j{E), i-e. using vector notation 

(B.29) r = R^e, with R = {rj^i)li^^,rj^i := ri{\j). 

The normalization condition for the reduces to X]f=i <^k,j'''j,i = ^k,i or equivalently 



(B.30) 



R = A-\ 



Note that the last three equations together with Lemma B.19 determine the holo- 
morphic differential completely. 

Finally the matrix of f3 periods of the ui 's can be expressed as 

g 



Tk,l 



T = BR = BA-^. 



Therefore 
(B.31) 

We combine equations ( ^^2^ ) , ( ^^251 ) , ( ^lE^ , ( ^^271 ) , ( |R3TD and Lemma pTg| to com- 
pute TTzr. 



diag (-2^(-1)s-'^)B(''^s) 



diag i27r{-iy-''hklnpk) + 

xC'diag(^(-ir'=). 
Introducing the matrices 

(B.32) n := diag Inpfc)^^^") diag ((-1)^"'=), 



(B.33) T2 := 



/ 1 

V 1 



diag (-7r(-l)e-i^)B('^^s)C diag ((-l)*^"'"), 



(B.34) 

we can summarize the information about vrir as follows. 
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Lemma B.35 Given the definitions made above, we can express 

irir = diag (Inp^) + t^'^^s) ^ 

where the entries of the matrix t^'^^s) (7I functions of the parameters (a, b, Xj ,Pj)- 
Furthermore they have an even extension for nonpositive values of pj , which pre- 
serves the regularity. 

Proof : It is only the matrix ri, which needs some consideration, namely the 
questions of a smooth and even extension for the terms {In pf:)cj^j\ But this is an 



immediate consequence of Lemma B.IS . 



O 



B.2 Frequencies and phase - U,V and Z 

As already remarked in Appendix ^ ( |A.l^ ), the vector U can be expressed as an 
integral of the uJi 's. More precisely 

u = — I J» 



1 




27ri J 


I 










J PS. 






2 



h VR(E) 

The last integration is performed on the upper sheet along the real axis. Using 
equations (|R2^ , and ( |R30D we conclude the following. 



Lemma B.36 Given the notation as above. Then 

(B.37) C/= C-i-,r^)C- (2 f ^.e). 

U is a smooth function of all the parameters {a,b,Xj,pj). Again there is an exten- 
sion for nonpositive values of the pj 's, which is even in each pj and preserves the 
smoothness. 
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Proof : The methods are the same as in all the other lemmas. Looking at 
the formula one might expect some difficulties for the differentiation with respect 
to b, but the transformation E ^ E -\- b resolves the situation. Secondly one might 
wonder whether the infinite domain of integration causes problems, but it is easy to 
check that all derivatives have uniformly at least as good a decay as which is 
integrable. 

O 



In order to evaluate V (see ( |A.8D ), we have to determine the zero order coefficient 
of at Poo and (see Q III.3.8 (3.8.2)]). Denote ri{E) = Y.'jZldijE^ and 
^ = as the local coordinates at the infinities. Then 

This implies 

uJi = {-di^g-i + O{0)dC at Poo, 
^/ = (t^i,3-i + 0(e))dCatPi. 



We obtain 



Introducing the abbreviation 

9 

(B.38) gk := J] 



m = 1 
m k 

it follows from (Q and (|B]2|) - ( |B3oD , that 

9 

i=i 

We can rewrite this equation as 
(B.39) V = -R^ diag {gt) 
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(B.40) 



y = -diag {-{-l)9-^)C^ diag {gu) 



I 1 



V 1 



The next lemma will finally tell us that we can choose the Aj 's as functions of 
the other parameters such that the solution of the Toda lattice which we have 
constructed in Appendix ^ is time periodic with frequency 7. 

Lemma B.41 There are smooth functions \j{a, b,pi, . . . ,Pg), 1 < J < 5, which are 
even in each pk such that 



(B.42) 



V 




Proof : This lemma is a consequence of the implicit function theorem. To see 
this more explicitly, we use Lemma B.l£ and equations ( |B.38 ) to evaluate 

the formula ( |B.40|) . 



V 



diag 



V 



m = 1 
m ^ k 



+ 0{p^) 



I 



(B.43) 

If all the Pk 's are equal to zero, we choose 



A,- 



a + b 



and remark below), in order to solve equation (B.42). Observe 



(A^''^ is real, see 

that the A^'''* 's are distinct and lie in (a, ^i^) C (a, 6). We compute the derivative 



5AF(p = 0,a,6,Af^) = diag 



-2Af ^ + a + & 



^.^[xf-a){h-Xf), 



which is invertible as k'^ < 2e~2 for all 1 < k < g and one can therefore choose 
ea and ej, in ( |B.2| ), (|B.3| ) small enough such that kj < ^ for all I < k < g. The 
evenness of the Xj 's in the pk 's is a consequence of the evenness of all the terms in 
the equation and the uniqueness of the solution. 



O 
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Remark B.44 We have chosen for each A^-°^ only one of the two possible solutions 



of the equation ( B.4^ ) for p = 0. We have seen in Remark 5. 3(\ , that the physical 



reason for this lies in the direction in which energy is transported in the correspond- 
ing g-gap solution. The reader may recall that exactly the same situation occurred 
in Chapter^ with the choice of the spatial frequency (3. 

Next we compute the phase Z. Therefore it is necessary to determine the vector 



of Riemann constants K. Proceeding as in |FK, VII. 1.2] we obtain 



9 

m=l 

as it is easy to check that the zeros of P ^ 'd{A{P)) are given by E2m~i, 1 < < 5- 
This imphes 

9 

(B.45) z = + 

m=l 




The path of integration from E2m-i to Pm is chosen along the real axis, beginning 
on the upper sheet and in case that Pm lies on the lower sheet we switch sheets at 
E2m- This choice is consistent with the description at the beginning of the proof of 



Theorem A.ll. 



The following lemma states that we can choose the Pj 's such that all phases Z 
are obtained. Recall that Pj is a point on Rg such that vr(Pj) G [E2j~i,E2j\, (see 
( B.7D ) i.e. each Pj lies on a cycle diffeomorhic to S"^. 



Lemma B.47 The map (Pi, . . . , Pg) i — > Z is a surjective map from {S^)^ to R^/Z^ 
for all choices of parameters {a,b, Xj,pj). 



Proof : By inspection of equation ( B.46 ) it is clear that Z G M (see Lemma B.19 ). 



It suffices to show that the image of the map is open and closed. First of all we 
have to convince ourselves that the map is differentiable. The only problem may 
occur at points Pm = E2m-i as the path of integration changes discontinously at 
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this point. More precisely, we have to investigate what happens if we add to the 



path of integration a cycle Um- By equation ( B.45 ) and the definition of oo we see 
that this adds a vector to Z, which is zero in R^/Z^. This settles the question of 
differentiability. Furthermore we know (see [FK, III. 11. 11 (Remark 1)] or [Gl, Thm 
17.20]) that the differential of the map has maximal rank g at any point. This shows 
that the image is open. That the image is closed is a consequence of the well known 
fact that the image of a compact set under a continuous map is again compact and 
hence closed. 



O 



Remark B.48 The zeros of S{A{-) +Un + Vt- Z\t). 

We now prove the asssertion, that the function P ^ ■d[A{P) + Un + Vt — Z\t) has 
exactly g one zero in each gap for all n E R, t G M , which was used in the proof of 
Theorem [A.ll| . By [ |FK| , VI.3.3, Theorem b] it suffices to show that there exists a 
choice of Pj with vr(Pj) G [E2j~i, E2j], such that 

Un + Vt- Z{Pi,...,Pg) = -Z{P[,...,P'g) in MP /Z^ . 

As [/, Z{Pi, . . . , Pg) are easily seen to be real valued vectors, the existence of such 



points follows from the Lemma |B.47| above. 

B.3 and 



lb 



Recall the definition of /g in ( A.15| ). Let 



This differential has the desired behavior at the poles and hence we only have to 
normalize in order to obtain 



(B.50) a;(i) = - ^ 
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Note that lo^^^ just changes sign if we switch from the upper to the lower sheet, and 
therefore Vs > : = /~ (see equation ( |A.15D ). We can express 



lim - ( u^^^ +lnE] 



^ (1) , 1. p\ , / /■ ..,(1) 1- 



The above integration takes place on the upper sheet. We recall, that the 5-gap 
solution should be brought into the form = cn+ small. Comparing with 

(|Alg| ) and ( [A3oD , we see that this implies a condition on the following quantity. 



(B.51) I:=/o+ + /o-. 

Lemma B.52 



(B.53) I = Jim -2 \^ ^0^ + InSj + < J Jq\u > 

is a smooth function of all the parameters (a, b, Xj,Pj) and has an even and smooth 
extension for nonpositive values of the pj 's. Furthermore there exists a unique 
smooth function b{a,pi, . . . ,Pg), which is even in each pk, such that 



(B.54) /(a,6,A„p,: 



c. 



It is understood that A,- is a function ofa,b,pk as determined in Lemma B.4I 



Proof : We break I into several pieces which we examine seperately. During 
this proof it is understood that all the square roots that appear take on positive 
values. Let us first examine the dependence of I on the parameters. 



-L 



V(A - o)(A - b) 



arccosh ( — ^ (e - 1 I + In 



b-a\ 2E J J \2 y4 4:{a + b - 2Ey 



which tends to In as — > cxd. 
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The remainder of the first term of / in equation ( B.53 ) is given by 



^{\-a){\-h) 



n 



1 



dX. 



(A-A,)^ / 



In order to faciliate the differentiaton with respect to 6, we shift the variable 
of integration by b. Furthermore it is not too comphcated to see that all the 
derivatives of the integrands with respect to the parameters are uniformly 
bounded by O(^) on the domain of integration and hence integrable. 



For the last term < "^o^^ U >, it is enough to cite Lemma B.3(; and to refer 



to the techniques introduced in the proof of Lemma B.l£. 



Our second goal is to determine h from the implicit function theorem. Let Aj depend 



on a,b,pk as given in Lemma B.41 . Then we can write / = I{a,b,pj), a smooth 
function satisfying 



-2 In. 



Let 
(B.55) 

then 



l{a,b,0) 

6^°) := a + 4e-t, 



d 



c. 

62 



db ^ ' ' 2 

The implicit function theorem yields the remaining claims of the Lemma. 



O 



Recall the definition of Rq in (A. 16). We proceed as above. Define 



(B.56) 

Then 
(B.57) 



(2) ._ 



(i^-^)nf=i(i^-A,) 



2/R(E) 



dE. 



UJ 



(2) 



UJ, 



(2) 



E 
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Equation ( |A.16[ ) yields 



— 



lim 

E~*oo 



lim 



-2 / ^(2) _ E 
Jb J 



-2 / LO, 

lb 



(2) 




E] + < I Jq\u> . 



The following lemma states that we can choose the parameter a in such a way, that 
the average speed of the solution will be 0. 



Lemma B.58 Denote 
(B.59) 

Then 



R :— Rq + Rq 



R = lim 



(2) 



2S I + 2 < / 4^\?7 > 



(B.60) 

E^oo \ Jb I J a 

There exists a unique smooth function a(pi, . . . ,Pg), which is even in each pk, such 
that 

(B.61) R{a,b,Xj,pj) = 0, 

where it is understood, that the Xj 's depend on a, b and the pk 's as it is determined 



in Lemma \B.4ii and b depends on a,pk as described in the previous Lemma. 



Proof : We proceed as in the proof of Lemma [B.52| . The only difference is 
that the first term in equation ( [B.60| ) is split up in 

rE 



b V(A-«)(A-6) 



-.dX - 2E 



and a remainder 



\ a+b 



^{X-a){X-b) 



n 



^=1 ^ 1 



dX. 



The first part tends to —{a + b), as E 00, whereas the second part can be 
dealt with as in the previous lemma. Furthermore this part vanishes identically for 
pi = ■ ■ ■ = Pg = 0. With the techniques of the proof of Lemma B.19 we observe 
that /q,^ Wq^"* = 0{p1). Finally we examine i? as a function of a and p^. From what 
was just said it is immediate, that 



R{a,pj = 0) = -(a + 6(a,0)) = -(2a + 46"^). 
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defining a^''^ := — 2e 2 we obtain 



i?(a(0),0) 



d 



— i?(a(0),0) = -2, 
which allows the use of the implicit function theorem. 



O 



Note that ( |A.19D , (|A.20D together with Lemmas |B.35| , |B.36|, |B.4lL |B.47| , |B.52| , 



B.58 already prove Theorem 5.8 with the exception of the last statement, concerning 



the limits pj ^ 0. 



B.4 The closing of gaps 

All the formulae of the last sections were derived under the assumption that all the 
the gaps are open, i.e. for all 1 < j < g : pj > 0. Nevertheless we have made it 
a point that almost all of the analytical expressions which we have calculated have 
a smooth (at least C^) continuation for nonpositive values of the pj 's. It is not 
a priori clear that the limits we obtain if we let some of the pj 's tend to zero, 
coincide with the formulae for the corresponding lower gap solution. It is the goal of 
this section to convince ourselves that the expression for ^(^(n,t) which is given in 



equation ( A. 19 ) has a continuous limit if some or all gaps close as described above. 
Clearly it suffices to examine the case that only one of the gaps closes at a time. 

To be more specific, fix pj > for 1 < j < g,j ^ i' and pi, = 0, as well as the 
other parameters a, b, \j, Pj. We further assume that g > 1. The case that the last 
gap closes will be dealt with at the end of this subsection. 

We use the following notation. All quantities which have to be evaluated will 
in two versions, namely with or without ' (e.g. A, A'). Without ' denotes the 
quantity for the g-gap case, in the limit pu — > 0, whereas the quantity with ' stands 
for the corresponding g — 1 gap expression with the same choice of the remaining 
parameters. Note that the parameters Xu,Pu and do not appear in the g — 1 gap 
case. 

For a g X g matrix M, we define by M^jj^u the (g — I) x (g — 1) matrix which is 
obtained from M by cancelling the -th row and column. Similarily for a vector v, 
we denote by Vj^i, the vector where the -th entry is cancelled. Finally the {k,j) 
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entry of a matrix M will sometimes be denoted by M{k,j). Let us now collect all 
the technical information, we will need. 

Proposition B.62 With the notation introduced above, the following relations hold 
for g > 1. 

(i) Akj^y = diag ( sgn {Xj, - Xu))A' diag {jj^)- 
M^J) = 0, for j / V. 

(a) Ckj^u = diag (Afc - Xu)C' diag ( sgn (Xk - Xy)). 
C{v,j) = 0, for j / u. 

C^^;^, = diag {Xk - A,)C"(H diag ( sgn (Afc - A,)). 
C('^")(z.,i) = 0, forj^v. 

(ill) Uj^, = U'. 

(tv) Vj^, = v. 

(v) Zj^y = Z'. 

(vii) I = I', 
(via) R = R'. 

Proof : 

(i) 



k = v : dkj 

(ii) We obtain the information about C from evaluating the relation AC = / in the 
following order. From the (z^, i^) entry conclude that C{i',i') = Looking 
at the v -th row we see that 'ij ^ v : C{y,2) = 0. This shows that Ckj^v is 



-.dE 



Xj — Xp J X 



k~Pk 



vwrn 



sgn {E - Xu)dE. 



Syj-Khy (see equation (|B.18D 
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the inverse of Akj^u, and this yields the claim for C. To prove the claim for 
it is enough (see (|B.21[) ) to observe from (|B.13D that for k ^ u 



irhk 



m 



For j 7^ u, we check that rr ^0^dE 



It 



°° ^ dE. The claim 



then follows from equation (|B.37 ) and from property (ii) (esp. C^{k^v) = 
for k ^ v). 

(iv) The proof is similar to (iii). Here it suffices by equation ( |B.4C1| ) to show that 

1 









diag {gk) 


01 


= dia 









Afc — Ajy 



diag {g'k) 



(v) 



But this follows directly from ( B.38 ). 

Looking at formula (B.46) it is sufficient to verify that 



-.dE 




3+^ 



diag ( 



At — Ai, 



■.dE 



It is not difficult to check this relation, if one observes that for j ^ v : 
fp" -^^dE = in the limit of Vv ^ 0. 



(vi) 



We have to examine the behavior of t\ and of T2 (see equations ( [B.32| ),( p.33| )). 
The fact that (t\)^ -^^^ = (ri)' follows directly from (ii) above and ( [B.32| ), using 
the fact that iov k u : hk = \\^^\ \h'k- The proof that {T2)k = (''"2)' needs 
more consideration. First we remark, that because of C{v,j) = for j ^ v 
(see (ii)), it suffices to prove the following equality. 

/ 1 \ 

: ■.. diag (-27r(-l)s-'^)B('''=s) 

Vi ... J 

1 . 

diag (-27r(-l)g-i-k) (B(-g))' diag (-^^ 
^ ^ Ak - At 



111 



Write 5(^*^9) = B- Observe that by equation (^1^ 



diag (-27r(-l)g-'^)B('^''^s) 



diag {27T{-iy-''hklnpk 



and one checks easily that this equals 
1 

diag (-27r(-l)S-i-'^) (b^'^'^s) 



diag 



Ak — A. 



To complete the proof we must check the relation for B itself. It suffices to 
compute bkj for j 7^ v. It is enough to verify the following two relations. 

• For A; / and A; / 1/ + 1: hj = -^^^^^pj^K'j'^ 
where k' = k, for k < u and k' = k — 1, for k > u, 
and j' = j, for j < u and j' = j — 1, for j > i^. 

• -bu,j + K+ij = xpi^Kj'- 

The first relation follows directly from the definitions. The crucial point for 
proving the second relation is the following calculation. In order to investigate 
the limit Pi/ — > explicitly, let us now assume for a moment that > 0. Let 
e > be an arbitrary but small number and let pi, < e. Then 



ejiE) 



dE + 



1 



\j - At 



e'^{E){E-K) 



-.dE 



+ 



e'^{E){E-K) 



T 



dE 



dE 



Xj - \ 



ds 



+ 



e'j{Xu + s) 



Xj — Xi/ 



e ^\R'{Xu + s) 



zds 



ds 



as pu 0. 
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(vii), (viii) 

The proof for / and R follows trivially from the formulae which were produced 



in the proofs of Lemma B.52 and Lemma |B. 58 , properties (iii) and (iv) of this 
proposition and the observation, that J^^ ujq^^ — > 0, as ^ for m = 1, 2. 



O 



Remarks: 

(1) In the above proposition we have kept a,b,Xj fixed as pu — > 0, but this is 
clearly not necessary. Indeed, if a, b, \j are given continuous functions of pu then 
the obvious analog of the proposition holds true. For example in the formula (vii), 

I{a,b,Xj,pi,...,0,...,pg) = I'{a,b,Xj, Pj^u), 

simply replace a, b, Xj by their limiting values at pi, = 0. 

(2) During the proposition we assumed that g > 1, as most of the terms do not 
carry any meaning for the 0-gap case. Only the quantities / and R are also well 
defined for g = and without changing the proof of property (vii) and (viii), we 
see that the proposition also holds in the transition from the 1 -gap to the 0-gap 
situation. 

We finally turn to the question of basic interest, namely how the formula for the 
solutions Xn{t) in Theorem |5.8| (see equation ( |5.1ClD ) behaves, if we let some or all 
of the Pu tend to zero. Recall that we have determined the parameters a, b and Xj 
as functions of pi, . . . ,pg such that x„(t) is a periodic function in t, which is of the 
form 

Xn{t) = cn + In ■ 



^{^U - Z\T)^{{n - + tV- Z\t) 



i9(-i[/ - Z|r)i?((n + i)C/ + tV- Z\t) 
Let us now again investigate what happens if one of the gaps closes, i.e. pu — > 0. 
The first question is, whether the choice of the parameters a, b and Xj as functions 
of the pk 's have the proper limit as Pi, — > 0, i.e. whether for example 

lim a(^>i, . . . ,p,y, ...,pg)= a'{pi, . . . ,Pu-i,Pu+i, . . . ,Pg). 

But this can be seen from properties (iv),(vii) and (viii) of Proposition [B.62| , as the 
solution of the determining equations ( |B.42D , ( |B.54D and ( [B.61j ) iov a = a{pi, . . . ,pg),b 
b(j)i, . . . ,pg) and Xj = Xj{pi, . . . ,pg) are unique. We are now in the position to apply 
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the Proposition B.62 in order to investigate equation ( 5.1C| ). It suffices to look at 
the behavior of the theta functions. Using Lemma B.35 and Lemma B.41 we obtain 

^{{n - + tV - Z\t) 
= Y.Pi ■■■■■ Pa' f 2vri(< l,U>{n- ^)- <l,Z >)+ < I, t^^-^s^I > 



(B.63) exp 



-i < I, 



> 7t 



V 



We see, that in the hmit py ^ Q only those terms in the sum survive, for which 
lu = 0. Using in addition the relevant results from Proposition B.62| we see that the 
limit of the i^-gap theta function is the appropriate g — 1-gap theta function. The 
other three theta functions in equation ( [5.10| ) can be dealt with in exactly the same 
way. Hence we have completed the proof of Theorem 5.8. 
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Appendix C 

Numerical experiments 



C.l Figures of lattice motion 



Figures C.6- C.11 below display the motion of the first ten particles (xq - xg with 
the zeroth particle on top and the nineth particle on the bottom of each figure) of 
lattices, which are described by the following system of equations: 

(C.l) Xn = F{Xn-l - Xn) - F{Xn - Xn+l), U > I, 

with driver xq of the form 

(C.2) xo{t) =t + e{sm-ft + 0.5cos2jt) 

and initial values given at t = 

(C.3) x„(0) = x„(0) = 0, n>l. 

We remark that we have also made experiments with driving particles xo(i) = 
2at + h{'yt) and h being periodic functions different from type ( p. 2)) (e.g. h piecewise 
linear) and we have always obtained results similar to those described below. The 
choice of parameters £,j,F is made as follows. On each page there are four figures. 
They correspond to different force functions F: 

top left : F{x) = (Toda lattice) 

top right : F{x) = 2.25x (linear lattice) 

bottom left : F{x) = 1.7l{x + 0.2^^) 

2.53 

bottom right : F(x) = 

1 — OAx 
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The parameters for these four types of force functions were chosen such that in the 
case of e = 0, the system behaves subcritical (i.e. 0.5 < acrit(-^)) and the lattice 
comes to rest {xn{t) —> cn as t —> cxd). For better comparison of the different force 
laws, we ensured in addition that F'{—c) ~ 2.25 in all for cases. Therefore we expect 
from the linear calculations (cf equation ( |1.26| )) that the threshold values for the 
frequencies should be approximately the same in all four cases (at least for small 
e), namely 7^ ~ 2>/k. Hence the different values which we selected for the driving 
frequency, 7 = 3.1,2.1,1.2 satisfy 3.1 > 71 > 2.1 > 72 > 1.2 > 73. Finally all 
experiments were made for two different values of the driver' s amplitude, namely 
e = 0.2 (see Figures ICj, IC^, |a|) and e = 0.5 (see Figures jal [ClO|, [OTlD. 



C.2 Spectral densities 



We consider equations (CJ.), ( C.3 ) in the case of the Toda lattice {F{x) = e^). 
As described in the Introduction we observe numerically that the spectrum of the 
corresponding Lax operator obtains a band-gap structure as t ^ cxd (cf Figures 1.29 
- |L3|). In Chapter | we derived under various assumptions an integral equation 
for the time-asymptotic spectral density 

{eigenvalues of L{t) that are < A} 



(C.5) 



J(A) 



lim 

t— >oo 



(cf (|1.34 ), Ansatz 2.22 and Theorem 2.271 ), which we solved explicitly in Theorem 
|2.3^ , provided a discrete number of data is known, namely the number and endpoints 
of the bands. In order to test the results of Chapter ^ we proceed as follows. We 
compute J(A) for A < inf aess{L{t)) numerically by evaluating ( p.Sj ) for large times t. 
This computation also yields a good approximation of the position of the endpoints 
of the spectral bands. Finally we determine the "predicted spectral density" by 
solving the integral equation given by ( 2.39D -( p.41 ), using the numerically obtained 
knowledge about the endpoints of the bands. In Figures |Cl^ and IC.13I below we 
display the results of these experiments for two different drivers xq of type ( p.2D , 
representing the cases 71 > 7 > 72 and 72 > 7 > 73, i-e. the one-gap and the 
two-gap situation. Both figures contain 

top left : the lattice motion 
top right : the time evolution of the spectrum of the corresponding 
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Lax operator 

bottom left : the numerically computed spectral density at t = 200 
bottom right : the "predicted spectral density" 

One observes very good agreement of numerical and predicted spectral densities, 
which a posteriori justifies the assumptions introduced in Chapter 0. 
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F{x) = (Toda lattice) 



F{x) = 2.25x (linear lattice) 




F{x) = 1.7l{x + 0.2x3) 

f\AVAr\rVAAVAr\AA/] 
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Figure C.6: Motion of lattices (cf (C.l) - (C.3)) with e = 0.2,7 = 3.1 
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Figure C.7: Motion of lattices (cf (C.l) - (C.3)) with e = 0.2,7 = 2.1 
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Figure C.8: Motion of lattices (cf (C.l) - (C.3)) with e = 0.2,7 = 1.2 
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Figure C.9: Motion of lattices (cf (C.l) - (C.3)) with e = 0.5,7 = 3.1 
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F{x) = (Toda lattice) 



F{x) = 2.25x (linear lattice) 




Figure C.IO: Motion of lattices (cf (C.l) - (C.3)) with e = 0.5,7 = 2.1 
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F{x) = (Toda lattice) 



F{x) = 2.25a: (linear lattice) 





Figure C.ll: Motion of lattices (cf (C.l) - (C.3)) with e = 0.5,7 = 1-2 
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Figure C.12: Toda lattice with driver xo{t) =t + 0.2 sin 7* + 0.1 cos 7*, 7 = 1.8, i.e. 
71 > 7 > 72; see Section C.2 for detailed description. 
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Lattice motion The spectrum 




Numerically observed spectral density Predicted spectral density 




spectral parameter, lambda spectral parameter, lambda 



Figure C.13: Toda lattice with driver xo{t) =t + 0.2 sin 7* + 0.1 cos 7*, 7 = 1.1, i.e. 
72 > 7 > 73/ see Section C.2 for detailed description. 
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Guide to the Video 



Part I and II of the video presentation follow rather closely Chapter 1 /resp. 
Chapter 2 of the text. In the third part, however, the viewer might find it diffi- 
cult to locate the details of what is discussed in the main text. Therefore we have 
provided this guide. It gives precise references to the main text (bold face) and 
repeats the formulae which were written on the black board for simple identification 
(underlined) . 

Part III 

We recall the long time behavior of the driven lattice and its dependence on the 
driving frequency 7, which was described in detail in Part I. The goal of this pre- 
sentation is to explain our result on the construction of ^ - time periodic solutions 
of the semi-infinite lattice, (cf Chapters 3-5), which are in good agreement with 
the numerically observed time asymptotic motion of the lattice. More precisely we 
construct solutions of 

(1) Xn = F{Xn-l - Xn) - F{Xn - Xn+l), U > I, 



with 



(2) xo{t) :=eY,h 



i-ymt 



(cf (3.1), (3.2)). For e = the lattice at rest Xnjt) = cn solves (0), (|). We ask 
whether these solutions persist for e 7^ 0? Our results can be summarized as follows. 



Let 7fc := fy/F'{-c), G N, be the sequence of threshold values for the driving 
frequency. 

• In the case of arbitrary force functions F (cf Section 3.2.2, General As- 
sumptions) and for 7 > 71 or 71 > 7 > 72 there exist eo > such that for all 
< e < eo we can construct ^ - time periodic solutions of (||) , (^) , satisfying 
Xnit) - cn = 0(e) (cf Theorems 3.56 and 4.24). 

• In the case of the Toda lattice {F{x) = e^), for all 7 G R+ \ {7^ : /c e N} there 
exist eo > such that for all < e < eo we can construct ^ - time periodic 
solutions of (H), @, satisfying Xn{t) — cn = 0(e) (cf Theorem 5.13). 
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It is instructive to consider the linear case F{x) = x. By separation of variables 
we find that 



(3) Xn{t) = eJ2h 



n i-ymt 



mi 



solves (D, (I) , if 

(4) zl, + 5mZm + 1 = 0, with 5^ = -2 + {jmf. 

We notice that a mode z'^e^"'"^^ corresponds to a travelling wave in the case \5m\ < 2, 
( Zm = e^^*" for f3m G K). For \8m\ > 2 one can choose \zm\ < 1, i-e. the correspond- 
ing mode is decaying exponentially in n. (cf (1.25)). The solutions we will construct 
in Chapters 3-5 will be nonlinear versions of (^). 

Our first result in Chapter 3 (Lemma 3.31 and Theorem 3.38) will reduce 
the construction of the desired solutions of (||), (|2|) to showing that there exists a 
sufficiently large parameter family of travelling wave solutions of the doubly infinite 
lattice. More precisely, in the case that the driving frequency satisfies 7fc > 7 > 7fc+i, 
we have to construct a 2k - parameter family of fc-phase waves. This reduction is 
proved by separating the equations for the "travelling wave part" and the "exponen- 
tially decaying part" (cf Lemma 3.31) and then by solving for the exponentially 
decaying part (cf Theorem 3.38) (see also Section 3.1). 

In the case of the Toda lattice the multi-phase travelling waves can be con- 
structed for arbitrary number of phases in terms of theta functions, using the well 
known g-gap solutions. Due to limited time these solutions are not discussed any 
further in this presentation. See Chapter 5, Appendix A and Appendix B for 
more detail. 

We will spend the remaining time to show how one can obtain single - phase 
waves in the case of arbitrary force functions F and what the difficulties are in the 
construction of multi-phase waves (in particular two-phase waves), which we have 
not yet overcome. By a two-phase wave we denote a solution of the doubly infinite 
lattice of the form 

(5) Xnjt) = cn + X2{nPi + 7t, n/32 + 27^), 

where X2 is 27r-periodic in both arguments. Functions of this form are ^ - time 
periodic. The spatial frequencies /3i,/32 turn out to be approximately given by the 
dispersion relation for the linearized lattice, 

(6) 2cos/3, = 2-ig^. 



2 



Remarks: 

(1) This formula was given incorrectly in the video presentation. 

(2) As mentioned above, two-phase waves are considered in the case 72 > 7 > 73, 
which means that there exist real solutions (3^ of (^) only in the case |A;| < 2. 

We first dicuss the construction of single-phase waves (see Chapter 4). Expand 

(7) Xnit) = cn + xi{n(3 + jt) 
in a Fourier series 

(8) xi(^) = E r{m)e^"'<^. 



The resulting equation for the sequence r can be written in the form 
(9) A(/3)r + W{r) = 0, 

where A(/3) is a linear operator depending on /3 and W{r) denotes the nonlinear 
term. The operator A(/3) is diagonal and its entries are given by 

A(/3)(m,m) = 2 cos /3m - 2 ^7!^^ , 

F'(-c) 

2 13m (7771)^ 



(10) = ¥+F'(-c) 

(See Lemma 4.10 for a derivation of (P), (px|)). 

Remark: 

In the video presentation we have given a different formula for A(/3) which can 
be obtained by dividing (^) by — 4sin^ which then represents the equation for 
the sequence f(m) := r(m)(e~*'^'" ~ !)• Furthermore we have assumed F'{—c) = 1, 
which can always be achieved by rescaling time. Therefore we obtain 

(7771)^ 



^11) A(/3)(m,m) = l 



4 sm ^ 



In the video we proceed to give a brief description of the Lyapunov - Schmidt 
decomposition which is explained in the introduction of Chapter 4 and after 
the proof of Remark 4.12. The final result is that we construct a two-parameter 
family of single phase waves of the doubly infinite lattice. 
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We repeat this ansatz in the case of two-phase waves, 



(12) Xnit) = cn + X2{npi + ^t,nl32 + 2-/t), 



(13) X2{ei,e2) = 

mGZ2 



i<m,d> 



The equation for f(m) := r(m)(e *</3'"^> — 1) can be written as 

(14) A(/?)r + W{r) = 0, 
with 

(15) A(g)H.^) = i- ^>;;,y . 



Note that we have again set F'{—c) = 1 without loss of generahty. Proceed- 
ing analogously to the single-phase case, we define the projection P, (Pr)(m) := 



ljl^l_^]^j.r(m), {\m\ = \mi\ + |m2|), Q := I — P. Equation (14) is equivalent to the 
following two equations. Denote /i := Pr, ip := Qr. 

(16) PA{P)n + PWifi + ip) = 0. 

(17) QAil3)^ + QWin + ip) = 0. 

The first step in the Lyapunov- Schmidt decomposition is to satisfy the infinite di- 
mensional equation ( p^ ) by choosing fi £ Ran (P) as a function of (3 and ip E Ran 
(Q) for (/3, ip) in a neighborhood of (/5(°), 0). We choose such that QA(/5(°)) = 0, 
i.e. /3W = (/3;°\/3f ), with 

(18) 4sin2^ = 7^ 



(19) 4sin2^ = (27)^. 



(0) 



Note again that (18) and (19) have real solutions (3!f^ because of the condition 
72 > 7 > 73- 

In order to solve (|l^) by implicit function theory, we have to investigate the 
invertibility of A(/3) : Ran (P) Ran (P). In contrast to the construction of the 
single-phase waves it is here that we encounter a major obstacle. In fact, a small 
divisor problem occurs, as there exist sites m G Z for which A(/3)(m, m) is arbitrarily 
close to zero. We call such sites singular sites. 
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We will determine the location of the sing ular sites of A(/3(o)) (cf Fig. || below). 
First of all one realizes that 72 > 7 > 73 implies that m £ 1? can only be a singular 
site if mi + 2m2 = ±1, ±2. We find the following cases. 

• For m G Z, satisfying m(-2/3f ^ + 4°^) « (mod 27r): 

(±1 — 2m, m), (— 2m, m it 1) are singular sites. They form a cross. 



For meZ, satisfying m(-2/?|°^ + 4°^] 
(—1 — 2m, m) is a singular site. 

For m G Z, satisfying m(-2/?f ^ + /3^°^) 
(1 — 2m, m) is a singular site. 

For m G Z, satisfying m(-2/?f ^ + 
(— 2m,m — 1) is a singular site. 

For m G Z, satisfying m(-2/?{°^ + 4°^] 
(—2m, m + 1) is a singular site. 



2/3j°^ (mod 2tt): 



2p[^^ (mod 2^): 



2/3^°^ (mod 2^): 



2^°^ (mod 27r): 




Figure 20: The location of the sing ular sites of A(/3W) 

A guide of how to construct solutions of (^) in such a situation is the work of 
W. Craig and C. E. Wayne on the existence of small solutions of the nonlinear wave 



equation |CW], where a similar small divisor problem occured. Their construction 
proceeds via a Newton iteration scheme. The crucial point is to control the small 
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eigenvalues of a truncation of the linearized operator. One recalls that in a Newton 
scheme it is necessary to linearize around the last approximation ti„ in order to 
obtain the better approximation Un+i and therefore the linearized operator is in 
general not in diagonal form with respect to the standard basis. The main technical 



tool in | CW ] is a technique introduced by J. Frohlich and T. Spencer( [FS| ) for 
the construction of inverse operators which makes it possible to first examine one 
singular site at a time and then to control the interaction between them. 

If one tries to apply this procedure to our problem one recognizes that the 



singular sites are not as well seperated as in |CW]. This causes additional difficulties 
which we have not resolved yet. Recent work of Surace (|S1], |^^) contains a 



version of the Frohlich-Spencer technique which might be helpful in overcoming 
these difficulties. 
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